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The  research  in  this  dissertation  is  motivated  by  the  basic  question  "What  can  and 

/ 

cannot  be  accomplished  by  feedback  control?"  In  particular,  this  dissertation  shall’ 
address  three  basic  issues:  a)  Determining  which  types  of  controllers  are  optimal  for 
certain  classes  of  control  problems,  b)  Investigating  the  absolute  limitations  of  feedback 
control  for  multiobjective  problems  and  the  performance  tradeoffs  available  between  the 
various  objectives,  and  c)  Developing  efficient  methods  for  synthesizing  robustly 
stabilizing  controllers  for  families  of  plants  featuring  block-structured  uncertainty. 

With  regard  to  the  issues  identified  above,  the  principal  contributions  of  this 
dissertation  can  be  outlined'as  follows.  First,  it  is  shown  that  for  the  problem  of  robustly 
stabilizing  a  family  of  plants  featuring  dynamic  uncertainty,  linear  time-invariant 
controllers  perform  as  well  as  arbitrary  nonlinear  time-varying  controllers.  Second,  a 
new  controller  synthesis  procedure  called  residue  iteration  is  developed  for  synthesizing 
robustly  stabilizing  controllers  for  families  of  plants  featuring  block-structured 
uncertainty.  This  method  is  simpler  and  numerically  more  attractive  than  any 
previously  existing  technique.  Finally,  an  algorithm  is  presented  which  enables  one  to 
compute  an  absolute  upper  bound  on  the  performance  levels  attainable  in  multiobjective 
H^-optimization  problems. 
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CHAPTER  1 


INTRODUCTION 


1.1.  Motivation  and  Background 

The  main  thrust  of  control  theory  since  the  1%0’s  has  been  in  the  direction  of 
controlling  accurately  modeled,  highly  structured  systems.  However,  in  many  situations 
available  models  are  partial,  inaccurate  and  far  from  being  clearly  structured.  These 
modeling  uncertainties  arise  from  a  variety  of  sources  such  as  linearization  of  nonlinear 
dynamics,  drift  of  plant  parameters  during  operation,  errors  due  to  model  reduction 
techniques,  and  identification  errors.  A  method  of  dealing  with  these  modeling 
uncertainties  is  to  allow  the  true  plant  model  to  be  any  of  a  family  F  of  possible  plant 
models.  Any  controller  designed  must  then  achieve  the  desired  performance  specifications 
not  only  for  the  nominal  plant  model,  but  rather  for  the  entire  family,  i.e.,  the  controller 
must  be  robust  with  respect  to  the  modeling  uncertainty  described  by  the  family  F.  In 
classical  control  theory  the  concepts  of  gain  and  phase  margin  are  employed  to  study  the 
robust  stability  properties  of  feedback  systems  (see  Horowitz  [1963]),  but  these  provide 
no  systematic  design  methodology  for  the  synthesis  of  robust  feedback  controllers. 

For  single-input,  single-output  systems  the  relationship  between  a  system’s 
robustness  properties  and  the  return  difference  matrix  is  well-known  (see  Bode  [1945]  or 
Horowitz  [1963]).  A  central  focus  of  robust  control  theory  over  the  years  has  been  an 
attempt  to  generalize  these  results  to  multiinput,  multioutput  systems  (see  for  example 
the  works  of  Cruz  and  Perkins  [1964],  Rosenbrock  [1974],  MacFarlane  and  Postlethwaite 
[1977],  Safonov  and  Athans  [1977],  and  Doyle  [1979]).  Recently,  there  has  been  a  growing 
interest  in  examining  robust  control  problems  via  frequency  domain  methods.  The  anpeal 
of  this  approach  may  be  linked  to  several  factors: 
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a)  For  many  control  system  design  problems,  frequency  domain  plant  models  are  often 
easier  io  obtain  than  state  space  plant  models.  Therefore,  it  is  useful  to  have  an 
alternative  design  methodology  available  for  problems  where  state  space  plant 
models  are  difficult  to  obtain. 

b)  Many  existing  state  space  design  techniques  assume  perfect  knowledge  about  the 
plant  model  and  provide  no  systematic  or  satisfactory  procedure  for  the  treatment  of 
robustness  issues.  This  is  at  least  partially  due  to  the  fact  that  state  space  models  are 
not  well-suited  for  representing  certain  types  of  modeling  uncertainty  such  as 
dynamic  uncertainty  (see  Chapter  2).  On  the  other  hand,  in  the  frequency  domain 
approach  this  uncertainty  can  easily  be  represented  by  a  ball  of  plant  models 
surrounding  a  nominal  frequency  domain  plant  model.  1'hus,  a  treatment  of 
robustness  can  be  readily  incorporated  into  both  the  problem  formulation  and  the 
design  procedure. 

c)  The  popular  Wiener-Hopf  design  procedures  make  restrictive  assumptions  on  the 
allowable  external  disturbance  inputs.  Stochastic  disturbance  inputs  are  modeled  as 
random  processes  with  known  means  and  power  spectral  densities.  The  resulting 
Wiener-Hopf  design  is  optimal  for  this  single  disturbance  input.  By  contrast,  the 
uniformly  (or  H,^-)  optimal  control  approach  to  this  problem  minimizes  the  worst- 
case  effect  of  a  set  of  disturbance  inputs  and  may  therefore  be  more  amenable  to 
certain  design  problems. 

In  response  to  these  concerns,  several  new  approaches  for  studying  robustness  issues 
have  emerged,  particularly  in  the  context  of  multiinput,  mullioutput  systems.  These 
include  uniformly  (or  !!„,-)  optimal  control  (see  Zames  [1981],  Zames  and  Francis  [1983], 
Francis  and  Zames  [1984],  and  O’Young  and  Francis  [1987]),  robustness  and  stability 
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margin  studies  (see  Cruz,  Freudenberg,  and  Looze  [1981],  Doyle  and  Stein  [ 1 98 1  ]. 
Lehtomaki  et  al.  [1981],  Safonov  [1980]  and  the  references  cited  therein),  /r-synthesis 
methods  (see  Doyle  [1984],  Doyle,  Wall,  and  Stein  [1982]  and  Safonov  and  Doyle  [1984]), 
and  robust  adaptive  methods  (see  Kosut  and  Friedlander  [1985]  and  Rohrs  et  al.  [1982]). 

From  a  conceptual  standpoint  the  frequency  domain  approach  has  introduced  some 
extremely  powerful  tools  for  control  system  design.  However,  due  to  the  efficiency  of 
state  space  computational  techniques,  almost  all  of  the  necessary  computations  involved 
in  synthesizing  controllers  are  still  performed  using  state  space  techniques.  In  fact,  an 
important  focus  of  current  research  is  on  development  of  efficient,  reliable  state  space 
techniques  for  synthesizing  controllers  which  achieve  desired  frequency  domain 
performance  criteria.  This  problem  is  further  addressed  in  Chapter  5. 

One  primary  focus  of  this  thesis  is  on  the  so-called  robust  stabilization  problem.  For 
recent  research  on  this  problem,  see  Safonov  [1980],  Doyle  and  Stein  [1981],  Kimura 
[1984],  Vidyasagar  [1985],  and  the  references  cited  therein.  This  problem  is  of  central 
importance  in  robust  control  theory  and  may  be  formulated  as  follows:  Given  a  family  F 
of  plant  models  which  represents  the  set  of  possible  "true"  plant  models  (accounting  for 
modeling  errors),  hnd  (if  possible)  a  single  feedback  controller  K  which  stabilizes  all 
plants  in  the  family  F.  In  complete  generality,  this  problem  is  extremely  difficult  and  no 
solution  is  known.  However,  for  certain  special  cases  of  interest  a  solution  may  be 
determined. 

For  instance,  Kimura  [1984]  has  obtained  a  complete  solution  to  the  problem  of 
robustly  stabilizing  a  family  of  plants  described  by  dynamic,  additive  uncertainty  using 
linear  time-invariant  controllers.  Also,  Tanne.nbaum  [1980]  and  Khargonekar  and 
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Tannenbaum  [1985]  have  studied  the  problems  of  stabilizing  one  parameter  families  of 
plants,  including  gain  and  phase  margin  problems,  using  linear  time-invariant  controllers. 

In  response  to  realistic  modeling  problems  such  as  those  encountered  in  decentralized 
systems,  Doyle  [1982]  has  introduced  the  notion  of  additive  block-structured  modeling 
uncertainty.  This  leads  directly  to  the  problem  of  robust  stabilization  of  families  of 
plants  featuring  additive  block-structured  uncertainty  which  reduces  to  a  /^-synthesis 
optimization  problem.  Such  problems  may  be  solved  iteratively  through  a  method  known 
as  A-Q  iteration  (see  Doyle  [1983]). 

While  robust  stabilization  is  a  key  feedback  system  design  attribute,  any  realistic 
problem  formulation  must  incorporate  performance  objectives  such  as  command  tracking 
and  disturbance  rejection.  A  particular  formulation  of  this  necessity  is  the  uniformly  (or 
Hm-)  optimal  control  problem.  Here  the  objective  is  to  minimize  the  worst-case  effect  of  a 
set  of  frequency  shaped  disturbances  on  the  output.  This  problem  formulation  differs 
from  the  popular  Wiener-Hopf  approach  to  controller  design  because  it  optimizes  over  a 
set  of  possible  disturbance  inputs  whereas  the  Wiener-Hopf  approach  optimizes  for  a 
single  disturbance  input.  A  more  complete  discussion  of  this  problem  is  contained  in 
Section  3.4. 

Based  on  this  brief  review  of  background  material,  it  is  now  possible  to  describe  the 
objective  and  contribution  of  this  dissertation. 

1.2.  Objective  and  Contribution  of  This  Dissertation 

The  research  in  this  dissertation  is  motivated  by  the  basic  question  "What  can  and 
cannot  be  accomplished  by  feedback  control?"  In  particular,  this  uissertation  shall 


address  three  basic  issues:  a)  Determining  which  types  of  controllers  are  optimal  for 
certain  classes  of  control  problems,  b)  Investigating  the  absolute  limitations  of  feedback 
control  for  multiobjective  problems  and  the  performance  tradeoffs  available  between  the 
various  objectives,  and  c)  Developing  efficient  methods  for  synthesizing  robustly 
stabilizing  controllers  for  families  of  plants  featuring  block-structured  uncertainty. 

In  addressing  the  issues  identified  above,  the  principal  contributions  of  this 
dissertation  can  be  outlined  as  follows.  First,  it  is  shown  that  for  the  problem  of  robustly 
stabilizing  a  family  of  plants  featuring  dynamic  uncertainty,  linear  time-invariant 
controllers  perform  as  well  as  arbitrary  nonlinear  time- varying  controllers.  Second,  a 
new  controller  synthesis  procedure  called  residue  deration  is  developed  for  synthesizing 
robustly  stabilizing  controllers  for  families  of  plants  featuring  block-structured 
uncertainty.  This  method  is  simpler  and  numerically  more  attractive  than  any 
previously  existing  technique.  Finally,  an  algorithm  is  presented  which  enables  one  to 
compute  an  absolute  upper  bound  on  the  performance  levels  attainable  in  multiobjective 
H^-optimization  problems. 

The  first  issue  examined  at  length  is  whether  or  not  nonlinear  time-varying 
controllers  provide  any  advantage  over  linear  time-invariant  controllers  for  the  problem 
of  robustly  stabilizing  a  family  of  plants  featuring  dynamic  uncertainty.  It  is  known  that 
for  certain  robust  stabilization  problems,  such  as  gain  margin  problems,  nonlinear  time- 
varying  controllers  are  far  superior  to  linear  time-invariant  ones  (see  for  example 
Khargonekar  and  Tannenbaum  (1985]).  This  leads  one  to  believe  that  nonlinear  time- 
varying  controllers  may  also  be  useful  for  other  robust  stabilization  problems.  However, 
it  is  shown  that  for  the  problem  of  robustly  stabilizing  a  dynamically  uncertain  family  of 
plant  models  this  is  not  the  case. 
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In  general,  adaptive  controllers  are  nonlinear  and  time-varying,  thus,  the  main  result 


described  above  demonstrates  that  adaptive  controllers  arc  not  useful  for  robustly 
*  stabilizing  families  of  plants  exhibiting  dynamic  uncertainty.  Due  to  the  proliferation  of 


controller  synthesis  methods,  results  of  this  type  are  invaluable  in  ruling  out  various 


synthesis  methods  for  certain  classes  of  control  problems. 
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The  above  result  is  not  entirely  unexpected  for  the  following  intuitive  reason.  Due 
to  the  nature  of  dynamic  modeling  uncertainty  (i.e.,  possibly  unbounded  McMillan  degree 
of  the  true  plant  model)  it  is  impossible  to  further  identify  the  true  plant  via  input- 
output  experiments.  Hence,  a  nonlinear  time-varying  (or  adaptive)  control  strategy  will 
not  prove  advantageous. 


I 
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The  development  preceding  this  main  result  rests  upon  several  key  intermediate 
results.  The  most  significant  of  these  results  demonstrates  that  the  small-gain  theorem  is 
both  necessary  and  sufficient  for  a  certain  class  of  nonlinear  operators. 

Another  issue  examined  in  detail  is  the  development  of  efficient  synthesis  techniques 
for  finding  robustly  stabilizing  controllers  for  a  family  of  plants  featuring  block- 
structured  uncertainty.  At  present,  the  only  available  technique  for  solving  this  problem 
is  the  A-Q  iteration  method  developed  by  Doyle  [1983],  From  a  numerical  standpoint, 
this  approach  is  complex  because  it  requires  the  solution  of  an  infinite  dimensional 
optimization  problem.  Further,  this  approach  has  some  drawbacks  with  regard  to  the 
computational  difficulties  in  determining  A  and  the  subsequent  need  to  find  a  rational 
approximation  for  A. 

In  response  to  these  difficulties,  a  new  controller  synthesis  procedure  is  developed  in 
this  thesis.  This  procedure,  based  on  a  key  observation  due  to  Safonov  [1986],  requires 
only  the  solution  of  a  finite  dimensional  optimization  problem.  An  added  benefit  of  this 
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approach  is  that  it  yields  new  insight  into  the  robust  stabilization  problem.  In  particular, 
clean  sufficient  conditions  on  the  nominal  plant  P0  which  guarantee  infinite  additive 
stability  margin  are  developed. 

A  crucial  step  in  recovering  the  controller  designed  via  the  residue  iteration  technique 
is  the  determination  of  an  outer  interpolating  function.  A  new  method  of  finding  such 
functions  is  presented.  This  method  yields  far  lower-order  outer  functions  than  those 
obtained  via  the  methods  of  Youla  et  al.  [1974]  and  circumvents  the  rational 
approximation  problems  which  arise  in  the  methods  of  Ball  and  Helton  [1979]. 

The  previous  discussion  has  focused  upon  a  single  objective  H^-optimization 
problem,  the  robust  stabilization  of  a  family  of  plants.  However,  many  interesting 
control  problems  reduce  instead  to  multiobjective  H^-opiimization  problems.  An  example 
of  this  is  the  problem  of  robust  stabilization  of  a  family  of  plants  with  optimal 
performance  for  the  nominal  plant.  These  multiobjective  problems  are  extremely  difficult 
to  solve  analytically.  Therefore,  the  development  in  this  thesis  centers  on  the 
presentation  of  a  systematic  algorithm  to  establish  absolute  bounds  on  the  attainable 
performance  levels.  These  bounds  are  established  graphically  and  the  resulting  curves  are 
useful  for  demonstrating  or  comparing  the  tradeoffs  inherent  between  various  competing 
objectives. 

1.3.  Organization  of  This  Dissertation 

The  remainder  of  this  dissertation  is  organized  as  follows.  Chapter  2  describes  in 
detail  various  paradigms  for  treating  the  problem  of  modeling  uncertainty.  Following 
this  in  Chapter  3  a  review  of  key  background  material  which  forms  the  foundation  of 
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work  in  this  dissertation  is  presented.  The  relative  merits  of  nonlinear  time-varying 
feedback  as  compared  with  linear  time-invariant  feedback  in  the  context  of  certain  robust 
stabilization  problems  are  investigated  in  Chapter  4.  Following  this,  Chapter  5 
concentrates  on  the  problem  of  synthesizing  robustly  stabilizing  controllers  for  families 
of  plants  featuring  block-structured  uncertainty.  Finally,  Chapter  6  describes  a 
systematic  algorithm  for  establishing  absolute  bounds  on  the  attainable  performance 
levels  in  multiobjecti^e  H^-optimization  problems. 

The  main  results  and  implications  of  this  dissertation  as  well  as  some  possible 
directions  for  further  research  are  summarized  in  the  final  chapter. 

1.4.  Notation 

R  =  field  of  real  numbers 

R+  =  set  of  positive  real  numbers 

C  =  field  of  complex  numbers 

C*  =  field  of  complex  numbers  -  10} 

D=  the  open  unit  disk  =  |z€C  :  1  z  I  <  1} 

eJ°  =  the  unit  circle  -  (z€C  :  I  z  I  =  1 } 

H  =  open  right  half  plane  =  {s€C  :  Re  s  >  0} 

H  =  closed  right  half  plane  =  js€C  :  Re  s  ^  0} 

LTI  =  linear  time-invariant 


LTV  = 


linear  time- varying 


NLTV  = 


nonlinear  time- varying 


CHAPTER  2 


MODELING  UNCERTAINTY 

Since  models  are  inherently  inaccurate,  key  tasks  for  the  control  system  designer  are 
to  analyze  and  reduce  the  effects  of  modeling  uncertainty  on  the  overall  system 
performance.  A  reasonable  formulation  of  this  problem  is  to  represent  the  true  physical 
plant,  not  by  a  single  nominal  plant,  but  instead  by  a  family  F  of  possible  plant  models. 
This  representation  of  modeling  uncertainty  falls  into  three  broad  categories. 

(a)  Parametric  Uncertainty.  In  this  case  P  consists  of  a  continuously  parameterized 
family  of  LTI  plant  models  as 

F  =  { P(i( z ) :  a  €  Q  C  Rn}  (2.1) 

where  the  parameter  vector  a  takes  on  values  in  some  compact  set  fi.  Note  that  this 
description  of  modeling  uncertainty  reflects  considerable  a  priori  knowledge  about  the 
plant  to  be  controlled.  This  is  because  implicit  in  (2.1)  is  the  presumption  that  the 
designer  has  been  able  to  isolate  the  sources  of  modeling  errors  in  the  plant  model  to 
specific  elements  such  as  gains  or  pole/zero  locations  and  to  conservatively  estimate  their 
ranges  by  compact  sets.  In  many  cases  adaptive  controllers  are  best  suited  for  treating 
this  type  of  uncertainty  (see  for  example  Astrom  [  1 983],  Kumar  [1985]  and  the  references 
cited  therein). 

(b)  Dynamic  Uncertainty.  This  type  of  modeling  uncertainty  is  described  by  a 
frequency  weighted  ball  of  perturbations  around  a  nominal  plant  model.  This 
representation  reflects  considerably  less  a  priori  knowledge  about  the  physical  plant  (as 
compared  to  parametric  uncertainty)  and  is  characterized  by  a  complete  lack  of 
knowledge  of  the  McMillan  degree  of  the  true  plant.  A  specific  instance  of  dynamic 
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uncertainty  is  the  additive  unstructured  uncertainty  as  in  Doyle  and  Stein  [1981] 
described  by 

F8  =  |P0(z)+iW,(z):  II  A  II  <  8}  .  (2.2) 

Here,  P0  is  a  nominal  LTI  plant  model.  A  is  a  normalized  perturbation  and  W,(z)  is  a 
stable  and  stably  invertible  frequency -dependent  weighting  function  which  characterizes 
the  relative  magnitude  of  modeling  uncertainty  at  various  frequency  levels.  Typically, 
W,(z)  takes  on  larger  values  at  higher  frequencies  reflecting  the  inability  to  model 
accurately  at  these  frequencies.  In  addition,  it  is  required  that  all  plant  models  in  F6  have 
the  same  number  of  unstable  poles  (if  A  is  LTI)  or  that  A  be  a  stable  NLTV  perturbation. 


£ 


In  many  problems  of  practical  interest,  the  representation  (2.2)  of  modeling 
uncertainty  is  U>o  conservative.  Motivated  by  this,  Doyle  [1982]  has  introduced  the 
notion  of  block-structured  uncertainty  wherein  only  certain  subsystems  of  the  plant 
model  exhibit  significant  modeling  errors.  Through  judicious  block-diagram 
manipulations  it  can  be  shown  (see  Doyle  [1982])  that  many  arbitrary  uncertainty 
structures  can  be  transformed  to  one  with  additive  block-structured  uncertainty.  Thus, 
one  may  focus  on  additive  block-structured  uncertainty  to  be  represented  as 


IP  s 

^t 

A, 

P0<z)  + 

W(z)  :  n  Aj  II  ^  5 

An 

(2.3) 


V,  -WJ 


Here  the  jth  uncertainty  block  A,  is  of  size  rjXr  and  the  assumptions  of  (2.2)  apply. 

(c)  Mixed  Uncertainty.  This  representation  of  modeling  uncertainty  involves  a 
combination  of  both  parametric  and  dynamic  uncertainties.  A  paradigm  for  this 


description  is 


Fs  =  |Po(z)  +  AWft(z)  :  II  A  II  ^8.  a€  Q  Q  Rn]  (2.4) 

and  clearly  represents  the  most  realistic  situation.  In  this  case,  the  accurately  modeled 
subsystems  of  the  plant  model  arc  described  by  a  parametric  representation  and  the 
poorly  modeled  subsystems  of  the  plant  model  (e.g..  the  high  frequency  dynamics)  are 
described  by  dynamic  uncertainty.  Recent  research  articles  (see  for  example,  Rohrs  et  al. 
[1982])  on  the  effects  of  unmodeled  high  frequency  dynamics  on  adaptive  control 
algorithms  address  this  type  of  modeling  uncertainty. 


CHAPTER  3 


DEFINITIONS  AND  BACKGROUND  MATERIAL 

This  chapter  outlines  a  mathematical  framework  which  forms  the  basis  for  the 
frequency  domain  approach  to  control  system  analysis  and  design.  Many  fundamental 
concepts  of  frequency  domain  control  are  reviewed  and  these  form  the  building  blocks  for 
much  of  the  further  development  in  this  thesis.  While  only  discrete  time  versions  of 
these  ideas  are  described,  continuous  time  counterparts  of  these  notions  may  readily  be 
obtained. 

This  chapter  is  organised  as  follows.  In  Section  3.1  basic  input-output  theory  is 
discussed,  focusing  principally  on  the  operator-theoretic  representation  of  systems. 
Section  3.2  treats  stability  theory  and  presents  a  well-known  parameterization  of  all 
stabilizing  controllers,  along  with  the  basic  elements  of  stable  coprime  and  inner-outer 
factorization  theory.  Finally,  in  Sections  3.3  and  3.4  two  principal  problems  in  frequency 
domain  control  are  introduced:  the  robust  stabilization  problem  and  the  uniformly  (or 
optimal  control  problem.  These  problems  shall  occupy  much  of  the  attention  of  this 
dissertation. 


3.1.  Input-Output  Theory 

Throughout  this  dissertation,  the  following  notation  from  Willems  [1971]  shall  be 
employed.  Let  Z+  be  the  sei  of  positive  integers  and  define  the  usual  Hilbert  space 

30 

12  =  {u  =  (  •  ■  •  u_ , ,u0.u, ,... ) :  Uj  £  R  ,  £  uj2  <  oo) 

i=— oo 
oo 

equipped  with  the  usual  norm  II  ul,=(  £  u^)*.  As  is  well-known,  1,  admits  a  direct- 
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sum  decomposition  as  12  =  h,  +  h2  where 


h2  =  }u  =  (ujj.Uj.u,,...)  :  Uj  €  R  ,  £  uf  <  °°} 


i=0 


Let  bh-,  be  the  unit  ball  in  h2.  Define  a  projection  operator  ttn  acting  on  the  space  of  all 
two-sided  sequences  (not  necessarily  square-summable)  of  real  numbers  as 

wN  :  (  •  •  •  u.j.Uq.u,,...)  -*(•••  u.j.Uq.u, . uN_,,0....) . 

As  in  Willems  [1971]  let 

h1(  =  |u  =  (u0,u1,u2,...):  u  €  h2  for  all  N  €Z+}  , 

i.e..  h-,e  is  the  space  of  all  locally  h,  sequences.  Similarly  the  space  l2e  can  be  defined  to 
be  the  space  of  all  locally  1,  sequences  of  real  numbers. 

An  m -input,  p-output  nonlinear  time-varying  (NLTV)  operator  F  :  h2™  — *  12PC  will  be 
called  proper  (causal)  if  7tn(u)  =  rrN(v)  implies  that  7tnF(u)  =  77-NF(v)  for  all  u,v  in  h2™ 
and  all  N  in  Z.  Note  that  for  a  causal  operator  F,  Im(F)  Q  h2P.  Also,  F  will  be  called 
strictly  proper  if  7rN(u)  =  ;rN(v)  implies  that  77rN+,F(u)  =  7rN+1F(  v)  for  all  u.v  €  h2™  and 
all  N  in  Z.  The  norm  of  F  is  defined  as 


F  II  =  sup 


u€h 


2 

a*0 


II  F(u)  II  , 


H  u  I 


(3.1) 


and  F  is  said  to  be  h-, -stable  (henceforth,  stable)  if  II  F  II  <oo.  In  addition,  define  II  F  II  bh^ 
to  be  the  norm  of  F  restricted  to  inputs  u  in  bh2.  Finally,  the  incremental  norm  of  the 
operator  F  is  defined  as 


II  F(u+v)— F(  v)  II  , 


(3.2 


11  F  11  me  =  SUP  - 

u ,  v  €  h  ?m  »  V  II  2 

U^(l 

and  F  will  be  called  incrementally  stable  if  II  F  II  inc<oo. 

If  F  is  a  linear  time-invariant  (LTD  operator,  then  it  may  be  associated  with  a 
transfer  function  matrix  F(x).  In  this  case  F  is  stable  if  and  only  if  F(z)  belongs  to  H^, 
and  it  can  be  shown  that  (see  Francis  [  1 987] ): 

I!  F  II  =  IFlz)l.  =  ess  sup  I  F(ei0)  I  .  (3  3) 

0^6$ 2n 

A  transfer  function  matrix  F(z)  will  be  called  unimodular  if  and  only  if  FtzjSH^  and 
F  !(z)  exists  and  belongs  to  HM. 

The  orthogonal  projection  f rom  1  ,p  onto  h,p  will  be  denoted  tt,  and  from  1  ,p  onto  (h2P)' 
by  77,.  Then,  with  any  NLTV  operator  F  one  can  associate  its  Hankel  operator  T(F)  and  its 
Toeplitz  ofierator  0(F)  defined  as 

T(F) :  h2™  — *  (h,p)A  :  u  — ♦  j72Fu 
0(F)  :  h,™  — *  h2P  :  u  — .  77,Fu  . 

This  definition  is  consistent  with  the  standard  time  domain  interpretation  of  a  Hankel 
operator  as  mapping  future  inputs  to  past  outputs,  i.e.,  the  Hankel  operator  defined  from 
the  anticausal  projection  of  F.  The  Hankel  norm  of  a  NLTV  operator  F,  denoted  by 
II  F  II  H,  is  defined  as 

PFKh=  lir(F)ll  .  (3.4) 

If  the  operator  F  is  time-invariant  it  is  easy  to  see  that  its  Hankel  norm  may  be 
characterized  alternatively  as 


AV- 


( 3.5  J 


0  F  II  H  =  lim  II  ^nF(I— 7tn)  II  . 

N~»oo 

For  a  finite  dimensional  LTI  operator  F,  its  associated  Hankel  operator  T(F)  is  of 
finite  rank  and,  consequently,  has  finitely  many  singular  values.  These  are  called  the 
Hankel  singular  values  associated  with  F.  and  are  fundamental  input-output  invariants  of 
a  system.  They  capture  the  complexity  of  a  system  model  vis-a-vis  various  feedback 
problems,  particularly  in  the  context  of  model  approximation.  It  is  clear  from  (3.4)  that 
II  F II  H  is  simply  the  largest  Hankel  singular  value  of  F.  Given  a  purely  nnticausal 
transfer  function  matrix  F(z),  its  Hankel  singular  values  may  be  computed  as  follows 
(see  Moore  [1978]).  Let  £(A,B,C,D)  be  a  canonical  realization  of  F(z).  Solve  the 
Lyapunov  equations 

AMAT-M  =  BBT  (3.6a) 

ATWA-W=CTC  (3.6b) 

for  the  (.ontrollability  and  observability  grammians  M  and  W,  respectively.  Then  the 
Hankel  singular  values  of  F(z),  denoted  by  o^F),  are  given  by 

Oj(F)  =  k  'h  (MW) .  i  =  1 .2 . N 

where  N  is  the  McMillan  degree  of  F(z). 

The  central  notion  of  an  adjoint  operator  is  now  defined. 

Definition  (3.1).  Let  H(  and  H,  be  Hilbert  spaces.  Let  T:H.— *H,  be  a  linear  operator. 

t  • 

The  adjoint  of  T  written  T  is  the  (it  can  be  shown  that  T  always  exists  and  is  unique) 
operator  T  :  H,— .11,  defined  by 

<y  ,  Tu>„  =  <T  y,u>n  Vu  €  H,  ,  v  €  H,  . 
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It  T  is  a  linear  time-invariant  operator  with  transfer  function  matrix  T(z),  then  its 

adjoint  T*(/)  has  the  transfer  function  matrix  T 1  ( — ). 

/ 

As  is  well  known,  the  set  of  all  NLTV  causal  operators  forms  a  left-distributive 
algebra  with  respect  to  operator  addition  and  composition.  This  algebra,  however,  is  not 
right-distributive,  i.e..  A( B+C)  ^  AB  +  AC  in  general,  and  this  is  a  principal  technical 
obstruction  in  dealing  with  NLTV  operators. 


3.2  Stabilizability  and  Factorization  Theory 

Consider  the  feedback  system  shown  in  Figure  3.1  where  P  and  K  are  causal  NLTV 
operators. 


FP— 9292 

Figure  3.1.  A  standard  feedback  system. 


Following  Desoer  et  al.[  1980]  and  Willems  [1971  ]  one  has  the  following 
Definition  (3.2).  Consider  the  NLTV  operator 


ui 

el 

d>(  P.K )  : 

— 

U2 

°2 

I  he  feedback  syslem  Figure  3.1  is  said  to  be  weil-posed  if  </>  is  single- valued  and  causal 
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Further,  a  controller  K  (possibly  NLTV)  is  said  to  (incrementally)  stabilize  the  plant  P  if 
the  operator  d>(P.K)  is  (incrementally)  stable. 

Remark  (3.1).  For  a  linear  plant  and  controller  pair  the  operator  <KP,K)  can  be  expressed 
in  its  usual  2x2  block  matrix  form 

(I+KP)-1  — KO+PK)-1 

<2>(P,K)  = 

P(I+KP)-1  (I+PK)"1 

However,  due  to  the  lack  of  right-distributivity  of  nonlinear  operators,  if  P  or  K  is 
nonlinear  then  d>(P,K)  cannot  be  expressed  in  this  form.  In  this  case  it  is  easy  to  see  that 
if  one  of  the  operators 

<£,  =  (I+KP)-1  <bi  =  (I-P(-K))"1 


&2  =  P(I+KP)  1  <t>4  =  -KU-P(-K)  1 

is  unstable, then  <2>(P,K)  is  unstable. 

□ 

Let  P0(z)  be  any  proper  LT1  plant.  Then  P0  admits  LTI  stable  coprime  factorizations 
as 


P0  =  ND  1  =  6  1  N 
XN  +  YD  =  I  .  NX  +  DY  =  1 


(3.7) 


where  N,  N,  D,  D.  X,  X,  Y,  and  Y  are  all  stable,  proper  LTI  operators.  A  given  plant  P0 
admits  many  similar  factorizations.  However,  all  of  these  factorizations  are 
unimodularly  related.  More  precisely,  the  following  result  may  now  be  obtained. 

Lemma  (3.1).  Let  P0(z)  be  a  proper  LTI  plant  with  a  LTI  stable  coprime  factorization  as 


in  (3.7).  Further  suppose  P0  admits  an  alternate  (possibly  NLTV)  stable  coprime 
factorization  as 


XjN,  +  Y,D,  =  I 

po=Ni^r1  =  R."1N1 

N,X,  +  D,  Y,  =  I 

Then. 


N  X 

N,  X, 

c 

o 

D  y 

D>  v, 

T 

> 

o 

X  Y 

u_I  o 

X,  Yj 

N  D 

o 

< 

N,  D, 

where  U  and  V  are  unimodular.  In  other  words,  all  stable  coprime  factorizations  of  a 
given  LTI  plant  are  unimodularly  equivalent  to  any  LTI  stable  coprime  factorization. 
Proof.  Choose  U  =  D,  'd.  Then  D,U  =  DjID^'d)  =  D  and  N,U  =  N^D,  *D  =  N.  Thus.it 
need  only  be  shown  that  U  =  D,  'd  is  unimodular.  Starting  with  the  identity 
XjNj+YjDj  =  I.  it  readily  follows  that 

X,N,U  +  YjD,U  =  XjNjDj-1  D  +  Y,D  =  X,N  +  Y,D  =  U 

which  implies  that  Similarly,  U  1  =  XN,-rYD,  and  u  Therefore,  U  is 

indeed  unimodular.  Choosing  V  =  NNt  1  and  employing  an  identical  argument  complete 
the  proof. 

□ 


The  following  theorem  repr'scnts  a  very  special  case  of  a  general  resu.t  due  to 
Ananthram  and  Desoer  [1984].  Their  results  give  a  complete  global  parameterization  of 


all  NLTV  stabilizing  controllers  for  the  class  of  strongly  stabilizable  nonlinear  plants. 
Their  results  may  be  specialized  to  LTI  plants  to  yield  the  following: 

Theorem  (3.1).  Let  P0(z)  be  a  strictly  proper  linear  time-invariant  plant.  Then  the  set  of 
all  NLTV  controllers  that  stabilize  P0  is  completely  parameterized  by 

K  =  (X  +  DQXY-NQ)-'  (3.8) 

where  the  Youla  parameter  Q  is  allowed  to  be  any  possibly  NLTV  stable,  causal  operator. 
Proof.  From  Khargonekar  et  al.[l985]  it  follows  that  any  LTI  plant  P0  can  be  stabilized 
by  a  stable,  periodic  linear  controller  K0.  Then  P0  admits  the  following  stable  coprime 
factorization 

P0  =  ND-1  =  D"'  N 

with 

XN  +  YD  =  I  ,  NX  +  DY  =  I 

where 

N  =  P0(I  +  K0P0)_1  .  D  =  (I  +  K0P0)_1  .  X  =  K0  ,  Y  =  I 

N  =  (1  +  P0K0)_1  P0  ,  D  =  (I  +  P0K0r'  ,  X  =  K0  and  Y  =  I . 

Using  results  from  Ananthram  and  Desoer  [1984],  the  set  of  all  NLTV  controllers  K 
which  stabilize  P0  can  be  parameterized  as 

K  =  K0  +  Qd-PjQ)-1  where  P,  =  N  =  P0(I  +  K0P0r’  (3.9) 

and  Q  is  allowed  to  be  any  stable,  causal  NLTV  operator.  This  is  now  simplified  as 

K  =  K0  +  Qd-NQ)-1  =  [K0(I— NQ)  +  QKl-NQ)"1 


=  [K0  +  (I-K0N)Q](1-NQ)  1  =  (K0  +  DQXl-NQ) 


(3.10) 


=  (X  +  DOXY-NQ)  1  • 


Now  assume  that  P0  has  an  alternate  stable  coprime  factorization: 

P0  =  N1D1'1  =  6“'  N, 

with 

X,N,  +  Y,D,  =  I,  N,X,  +  D,Y,  =  1  .  (3.11) 

Fmploying  this  factorization  the  set  of  all  controllers  that  stabilize  P0  are  parameterized 
as 

K,  =(X, +D1QI)(Y,-N1Qir1  (3.12) 

where  Q,  is  allowed  to  be  any  stable,  causal  NLTV  operator.  However,  from  Lemma  (3.1 ) 
it  follows  that  all  coprime  factorizations  of  a  strictly  proper  LT1  plant  are  unimodularly 
related.  Thus,  without  loss  of  generality  X,  ,  D,  ,  Y,  ,  and  N,  may  be  chosen  from  any 
LTI  stable,  coprime  factorization  of  P0  completing  the  proof. 

□ 

Remark  (3.2).  The  set  of  all  LTI  controllers  K  that  stabilize  P0  may  be  parameterized  by 
(3.8)  by  restricting  Q  to  be  any  LTI  stable,  causal  operator.  Also,  observe  that  a 
controller  K  incrementally  stabilizes  P0  if  and  only  if  the  corresponding  Youla  parameter 
Q  is  incrementally  stable.  Further  from  (3.7)  and  (3.8)  it  follows  that 


(I  +  P0K)_I  =  (Y  —  NQ)  6 

(3.13a) 

K(1  +  P0Kf 1  =(X  +  DQ)  6 

(3.13b) 

P0K(I  +  P0K)"'  =N(X  +  QD). 

(3.13c) 

□ 

It  is  significant  to  note  that  all  three  of  these  transfer  functions  are  affine  in  Q,  i.e..  they 


all  have  the  general  form  A+BQC  where  A,  B,  and  C  are  fixed  stable  proper  transfer 
functions  (see  also (3.2 1 ) ). 

□ 

Consider  the  general  feedback  system  shown  in  Figure  3.2. 


Figure  3.2.  Feedback  system  for  the  small-gain  theorem. 

A  central  theorem  due  to  Sandberg  [1964]  and  Zames  [1966]  concerning  the  stability  of 
this  system  is  the  following 

Theorem  (3.2).  (The  Small-Gain  Theorem).  Consider  the  feedback  system  shown  in 
Figure  3.2  and  assume  that  M  and  A  are  causal  operators.  If  II  M  II  ^  I  then  the  feedback 
system  is  stable  for  all  operators  A  with  II  A  II  <1. 

Remark  (3.3).  It  is  crucial  to  note  that  the  small-gain  theorem  provides  only  a  sufficient 
condition  for  the  stability  of  a  feedback  system.  For  certain  special  cases  (i.e.,  when  M 
and  A  are  both  LTI  operators),  it  is  known  that  the  small-gain  theorem  is  both  a  necessary 
and  sufficient  condition  for  system  stability.  Indeed,  a  key  result  of  this  dissertation  is 
that  the  small-gain  theorem  is  both  necessary  and  sufficient  for  stabilizing  a  system 
consisting  of  a  particular  class  of  NLTV  operators. 


The  powerful  tools  of  inner-outer  factorizations  are  now  introduced.  For  a  more 
detailed  discussion  refer  to  Duren  [1970]  or  Francis  [1987],  A  matrix  transfer  function  F 
in  RH^  is  inner  if  FF  =  1.  Inner  transfer  functions  are  all-pass.  In  order  for  F  to  be 
inner,  clearly  F  must  be  tall  (number  of  rows  ^  number  of  columns).  An  important 
property  of  inner  functions  is  that  they  preserve  inner  products  under  left-multiplication, 
i.e.,  for  any  Ge  L^, 

I  FG  11=  II  Gil  . 

oo  oo 

A  matrix  transfer  function  F  in  RH^,  is  outer  if  for  every  z  with  Izl  <1,  F(z)  has 
full  row  rank,  i.e.,  F  has  a  right-inverse  which  is  analytic  in  Izl  <1.  Outer  functions 
must  be  wide  (number  of  columns  ^  number  of  rows). 

Every  matrix  transfer  function  F^RH^  admits  an  inner-outer  factorization  as 

F  =  F;F0  (3.14) 

where  Ft  is  inner  and  F0  is  outer.  These  factorizations  may  be  computed  using  spectral 
factorization  theory  and  involve  the  solution  of  a  set  of  algebraic  Riccati  equations  (see 
Doyle  [1983]  or  Francis  [1987]). 

In  addition,  a  matrix  transfer  function  F  in  RH^  is  said  to  be  co-inner  or  co-outer  if 
FT  is  inner  or  outer,  respectively.  Every  matrix  transfer  function  FCRH^  admits  a  co- 
outer-co-inner  factorization  as 

F  =  F„P„  (3.15) 

where  Fco  is  co-outer  and  FC1  is  co-inncr.  These  factorizations  arc  obtained  from  an  inner- 
outer  factorization  of  FT.  An  important  property  of  co-inner  functions  is  that  they 
preserve  inner  products  under  right-multiplication,  i.e.,  for  any  G€  L^, 
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II  GF  II  „  =  II  G  II  . 

3.3.  The  Robust  Stabilization  Problem 

As  described  earlier,  plant  modeling  uncertainty  may  be  treated  by  allowing  the  true 
plant  model  P  to  be  any  of  a  family  Fs  of  possible  plant  models.  Thus,  instead  of 
stabilizing  a  single  plant,  a  controller  must  stabilize  an  entire  family  of  plant  models. 
This  leads  naturally  to  the  robust  stabilization  problem  which  may  be  formulated  as 
follows:  Given  a  family  of  plants  F5  find  (if  possible)  a  single  controller  K  which 
stabilizes  every  plant  in  the  family  F5.  If  such  a  controller  exists,  then  it  is  called  a 
robustly  stabilizing  controller  and  the  family  of  plants  is  called  robustly  stabilizable. 

The  robust  stabilization  problem  has  been  extensively  studied  for  families  of  plants 
F5  featuring  unstructured  uncertainty  as  described  in  (2.2)  (see  for  example  Doyle  et  al. 
[1982],  Doyle  and  Stein  tl98l],  Glover  [1984],  Kimura  [1984],  and  Verma  [1985]).  Kimura 
[1984]  has  obtained  necessary  and  sufficient  conditions  for  robust  stabilizability  of  F5 
using  LTI  controllers  based  on  classical  Nevannlina-Pick  interpolation  theory.  Of 
particular  relevance  to  this  dissertation  is  the  following  result  due  to  Doyle  and  Stein 
[1981]  and  Chen  and  Desoer  [  1 982]. 

Theorem  (3.3).  Consider  a  family  of  plants  Fs  as  described  in  (2.2)  and  introduce  stable 
coprime  factorizations  of  P0  as  in  (3.7).  Then  the  following  are  equivalent: 

(a)  There  exists  a  LTI  controller  K  that  robustly  stabilizes  Fs. 

(b)  inf  II  W2K(I+P0Kr‘  <  -. 

LTI  K  stabilizing  P0  8 


CcJ  inf  II  w'jXD  +  W,DQD  II  ^  <  — . 

Q€H.  *  "  8 


The  corresponding  robustly  stabilizing  LT1  controllers  may  be  readily  computed 
using  a  variety  of  methods  (see  for  example  Vidyasagar  [1985]). 

Consider  again  the  family  of  plant  models  Fs  and  define  PU=(P0W,  ’)+  and 
Ps  =  P0W,  Pu-  Let  N  =  McMillan  degree  of  Pu.  Working  independently  Glover  [1986] 
and  Verma  [1985]  have  obtained  the  following  alternate  characterization  of  robust 
stabilizability  of  the  family  Fs: 

Theorem  (3.4):  There  exists  a  LTI  controller  K  that  robustly  stabilizes  the  family  Fs  if 
and  only  if 


o(r(P0W,  *))  >  8. 


(3.16) 


In  this  event,  a  robustly  stabilizing  K  may  be  determined  by 

K  =  -wr'(PN_,  +PS)-1 

where  PN_,  is  any  (N-l)st  order  optimal  HankeJ  norm  approximant  of  Pu  such  that 
Pu— PN_,  is  all-pass,  i.t.. 


Pu-Pn-,»« 


=  «npu-pN_1)ii  =a(r(p0w2 ')) . 


Remark  (3.4).  It  is  important  to  note  (see  Glover  [1986]  and  Verma  [1985])  that 
Theorem  (3.4)  holds  with  the  relaxed  assumption  that  all  plants  in  F5  have  only  the  same 
or  fewer  poles  as  the  nominal  plant  P().  However,  applying  some  basic  results  from  model 
approximation  theory  (see  Glover  [1986]).  it  readily  follows  that  even  under  this  relaxed 
assumption  the  family  Fs  is  robustly  stabiiizable  if  and  only  if  ail  plants  in  F5  have  the 
same  number  of  unstable  poles.  Therefore,  Theorems  (3.3)  and  (3.4)  are  indeed 


interchangeable  robust  stability  criteria  for  families  of  plants  Fs  as  described  by  (2.2). 


Motivated  by  these  results  and  by  more  general  synthesis  problems  such  as  robust 
performance,  Doyle  [1982]  has  introduced  motion  of  a  structured  singular  value  m 
which  provides  a  method  of  tackling  the  robust  stabilization  problem  for  families  of 
plants  featuring  block-structured  uncertainties.  The  theory  developed  in  Doyle  [1982]  is 
particularly  clean  for  the  case  of  three  or  fewer  blocks  (n^3)  where  an  exact 
determination  of  /j.  is  possible  via  the  A-scaling  iteration  (see  also  work  by  Fan  and  Tits 
[1986]).  The  following  development  concentrates  on  this  case  (n^3)  and  for  the  case  of 
more  than  three  blocks,  the  subsequent  results  yield  conservative  estimates  of  robust 
stability  margins. 

Following  Doyle  [1982]  define  the  set 


Cl  = 


0 


O 


d2hi 


dnIin 

n  jn 


(3.17) 


where  Ij  is  the  jxj  identity  matrix.  The  following  result,  analogous  to  Theorem  (3.3),  is 
easily  obtained. 

Theorem  (3.5).  Consider  the  family  of  plants  F5  as  described  in  (2.3)  with  block- 
structured  uncertainty.  Introduce  stably  coprime  factorizations  of  P0(z)  as  in  (3.7).  Then 
the  following  are  equivalent: 

(a)  There  exists  a  LTI  controller  K  that  robustly  stabilizes  Fs. 

(b)  inf  II  AW,K(I  4  PcKr'  A"'  11  M  . 

K  stabilizing  P0  8 

u  n 


(3.18) 


(c)  inf  II  A(W,XD  +  W,DQD)A  '  II  M 

Q€H„  '  "  °°  8 

V  f 

□ 

An  approach  lo  solving  the  optimization  problem  suggested  by  (3.18)  of  the  above 
theorem  is  the  A-Q  iteration  (or  /x-synthesis)  method  of  Doyle  [1983]  and  Doyle,  Wall, 
and  Stein  [1982].  Essentially  this  method  proposes  computing  the  infimum  in  (3.18) 
iteratively  by  freezing  A,  computing  the  optimal  Q,  freezing  Q,  and  finding  the  best  new 
scaling  A  and  so  on.  However,  although  the  optimization  problem  (3.18)  is  convex  in  both 
Q  and  A  it  is  not  jointly  convex  (see  Doyle  [1983]).  Therefore,  there  is  no  guarantee  that 
the  A-C  iteration  method  will  converge  to  a  global  minimum.  In  addition,  the  A-Q 
iteration  method  is  numerically  unappealing  due  to  the  complexity  of  computing  A  and 
then  having  to  find  a  suitable  rational  approximant  for  A  as  required  during  each  step  of 
the  iteration. 

In  Chapter  5  a  key  observation  due  to  Safonov  [1986]  is  employed  to  develop  an 
alternative  to  A-Q  iteration  for  finding  robustly  stabilizing  controllers  for  families  of 
plants  featuring  block  structured  uncertainty.  This  alternative  procedure,  called  residue 
iteration,  is  numerically  attractive  in  comparison  with  A-Q  iteration  because  it  replaces  an 
infinite  dimensional  optimization  problem  with  a  finite  dimensional  optimization  problem. 
Furthermore,  residue  iteration  is  relatively  easy  to  implement  numerically,  and 
circumvents  the  numerical  difficulties  involved  with  computing  A. 
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3.4.  Uniformly  or  H^-Optimal  Control 

In  addition  to  feedback  system  stability,  in  practical  control  problems  it  is 
imperative  to  incorporate  some  performance  objectives  such  as  disturbance  rejection  or 
reference  command  tracking.  1'his  objective  falls  naturally  within  the  framework  of  the 
following  uniformly  (or  H^-)  optimal  control  problem.  Consider  the  system  shown  in 
Figure  3.3. 


FP— 9290 

Figure  3.3.  Feedback  system  for  the  uniformly  optimal  control  problem. 

Here,  P0  represents  a  LTI  plant  model,  d  is  a  disturbance  input  of  unit  energy  and  W,  is  a 
coloring  filter  (weighting  function)  which  shapes  the  frequency  characteristics  of  the 
disturbance  input.  Without  loss  of  generality  Wj  can  be  chosen  to  be  stable  and 
minimum  phase.  The  objective  in  this  problem  is  to  design  a  controller  K  which  stabilizes 
P0  and  minimizes  the  worst-case  effect  of  any  disturbance  d  of  unit  energy  on  the  output 

y- 

By  straightforward  calculation  it  follows  thal 


y  =  (1  +  P0K)  W,d 


(3.19) 


and  hence  the  objective  of  the  uniformly  optimal  control  problem  is  to  minimize 
11(1  +  P0K)  'w,  II  oo.  From  (3. 13a)  this  problem  readily  reduces  to 


"J*  T  &  *>  "m 


/A 
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inf  MY— NQ)DW,  IM  (3.20) 

Q€RH„ 

and  the  optimal  controller  is  obtained  by  substituting  the  minimizing  Q  in  (3.20)  into 
(3.8). 

In  the  special  case  of  LTI  controllers, both  the  robust  stabilization  problem  described 
in  the  previous  section  and  the  uniformly  optimal  control  problem  described  above  reduce 
to  the  same  optimization  problem 

inf  II  A  +  BQC  II  „  (321) 

Qsrho> 

where  A,  B,  and  C  are  fixed  stable  proper  transfer  functions.  To  further  interpret  this 
problem  consider  the  system  shown  in  Figure  3.4. 


F  P—9288 


Figure  3.4.  The  model-matching  problem. 

Essentially  the  objective  of  the  optimization  problem  (3.21)  is  to  select  Q  in  HM  to  fit 
A+BQC  into  a  disk  of  smallest  radius  in  HM.  Hence,  (3.21)  is  sometimes  referred  to  as  a 
one-disk  problem.  Alternatively,  (3.21)  may  be  considered  as  a  proolcm  of  choosing  the 
stable,  proper  transfer  function  Q  such  that  the  transfer  function  -BQC  provides  the  best 
L^-norm  approximate  of  the  transfer  function  A.  Hence,  this  problem  is  also  referred  to 
as  a  model-matching  problem.  Any  solution  to  (3.21)  possesses  the  following  "all-pass" 


property  (see  Zames  and  Francis  [1983]). 

Theorem  (3.6).  Let  Q*  be  a  solution  to  the  one-disk  problem 

min  I!  A  +  BQC  II  ^ 

Q€RHo> 

where  A,  B.  and  C  are  RH^- functions.  Let  y  =  II  A  +  BQ*C  II  eo.  Then  for  any  frequency 
<a€R,  II  (A  +  BO*C)(jo>)  II  =  y. 

The  one-disk  (model-matching)  problem  can  easily  be  converted  into  the  well-known 
General  Distance  Problem  as  follows  (see  Francis  [1987]).  Consider  again  the 
optimization  problem  (3.21 )  and  let  B  =  BjB0  and  C  =  CcoCci  be  inner-outer  and  co-outer- 
co-inner  factorizations  of  B  and  C.  respectively.  Define  the  operators  E,  and  E2  as 

E,  =  and  E,  =  (Cc*  I-Cc'Cci) .  (3.22) 

l-B.B,* 


It  is  easy  to  show  that  E,  is  inner  and  E2  is  co-inner.  These  operators  are  called  dilations 
(see  Doyle  [1984])  and  possess  special  norm-preserving  properties.  Using  (3.22)  the 
optimization  problem  (3.21)  can  be  rewritten  as 


inf  II  A  +  BQC  II  =  inf  II  E,  ( A  +  B^OC^C^)  E2  II 
Q€RH„  QeRH„ 


=  inf 

T„-Q 

t21 

9€RH„, 

ti2 

t22 

(3.23) 


where  T„  =  B,  ACcj  .  T12  =  (I-B,Bi  )ACc‘.  T,,  =  B,  A(I-CCICcj).  T22  =  (I-BiBi  )A(I-CclCc|) 
and  Q  =  B0QCco-  This  problem  is  known  as  the  4-block  problem. 


The  solution  to  the  4-block  problem  is  iterative  in  nature  and  is  commonly  referred 
*o  as  y-iteration  (sec  for  example  Chu  [1985],  Chu  and  Doyle  [1985],  and  Francis  [1987]). 


In  certain  cases  the  dilation  matrices  E,  and  E,  introduced  in  (3.22)  are  trivial.  If  B  is  a 
square  full-rank  matrix  then  Bi  B,  =  BB,*  =  I  and  hence  E,  =  ( B,  O)*.  Similarly,  if  Cc  is 
square  and  t  ull  rank,  then  CflCcl  =  CclCcj  =  I  and  hence  E,  =  lCcl  O).  Whenever  one  of  the 
operators  B,  and  Cci  are  square  and  full  rank  then  the  4-block  problem  of  (3.23)  may  be 
reduced  to  a  2-block  problem  of  the  form 


or 


inf 

Q€RHo, 


t„-Q 

T„ 


(3.24) 


inf  llTn-Q  T,,  II  _  . 

Q€RHo 

These  problems  arc  also  solved  via  the  methods  of  y-iteration. 

If  however,  both  of  the  operators  Bt  and  Cci  are  square  full-rank,  then  the  4-block 
problem  of  (3.23)  reduces  to  a  basic  approximation  problem  of  the  form 


inf  II  G  —  H  II 

HCRH, 


(3.25) 


•  * 


where  G  =  B(  ACc,  and  H  =  B0QCco.  This  problem  can  be  interpreted  as  one  of  Ending  a 
purely  stable  transfer  function  matrix  H  which  best  approximates  (in  L^-norm)  the 
unstable  transfer  function  matrix  G.  It  should  be  noted  that  for  the  problem  of  uniformly 
(or  Hw-)  optimal  control  of  "short-fat  plants,"  i.e.,  plants  with  more  inputs  than  outputs, 
or  for  the  problem  of  robust  stabilization  of  families  of  plants  featuring  additive  or 
multiplicative  representations  of  dynamic  uncertainty,  the  matrices  B;  and  Cci  will 
always  be  square. 


The  basic  La,  H„  approximation  problem  stated  in  (3.25)  is  deeply  connected  to 


inU  '-polauon  theory  and  the  classical  moment  problem  and  has  been  studied  by  several 
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researchers  (see  for  example  Krein  and  Nudeiman  [1977],  Nehari  [1957],  Nevanlinna 
[1919],  and  Sarason  [1907]).  A  particularly  elegant  result  due  to  Adamjan  et  al.[1978]  is 
the  following: 

Theorem  (3.7).  Consider  the  optimization  problem  of  finding  a  closest  H^-matrix  H  to  a 
given  L^-matrix  G  as  stated  in  (3.25).  There  exists  a  solution  to  this  problem  and  this 
solution  is  characterized  by 

IIG-H  II  „  =  II  T(G)  II  .  (3.26) 

□ 

Remark  (3.5).  In  general  the  optimal  H  which  satisfies  (3.26)  is  nonunique.  In  basic 
terms  Theorem  (3.7)  states  that  the  distance  between  a  given  noncausal  system  G  and  the 
nearest  causal  system  H  is  given  by  the  norm  of  the  Hankel  operator  HG)  associated  with 
G.  This  theorem,  however,  does  not  provide  any  insight  into  determining  the  optimal 
H^-matrix  H.  Methods  for  finding  the  optimal  H^-matrix  H  have  been  developed  by 
Kimura  [1984],  through  classical  Nevannlina-Pick  interpolation  theory,  and  by  Glover 
[1984]  via  the  theory  of  Hankel  operators. 

□ 


CHAPTER  4 


NONUNEAR  TIME-VARYING  CONTROLLERS  FOR  ROBUST  STABILIZATION 

I  his  chapter  is  concerned  with  the  problem  of  robust  stabilization  of  a  family  of 
plants  exhibiting  dynamic  uncertainty.  Knnura  [1984]  has  obtained  a  complete  solution 
to  this  problem  in  the  case  of  unstructured  additive  uncertainty,  in  the  context  of  linear 
time-invariant  feedback.  The  intent  of  this  chapter  is  to  investigate  the  relative  merits  of 
nonlinear  time-varying  controllers  for  robustly  stabilizing  a  family  of  plants  of  the  type 
treated  by  Kimura  [1984].  It  has  been  demonstrated  in  Khargonekar  et  al.  [  1 985]  that  for 
certain  robust  stabilization  problems,  such  as  gain  margin  problems,  nonlinear  time- 
varying  controllers  are  far  superior  to  linear  time-invariant  ones.  This  result  suggests  that 
NLTV  control  may  ofler  advantages  in  other  robust  stabilization  problems  as  well. 
However,  it  is  demonstrated  in  this  chapter  that  for  the  problem  of  robust  stabilization  of 
families  of  plants  with  dynamic  uncertainty,  this  is  not  the  case.  Intermediate  to  proving 
the  above  theorem  are  several  results  themselves  of  interest.  In  particular,  it  is  shown 
that  the  small-gain  theorem  is  both  necessary  and  sufficient  for  a  certain  class  of 
nonlinear  operators. 

The  results  of  this  chapter  are  presented  for  discrete-time  systems.  However,  all  of 
the  results  very  easily  extend  to  continuous-time  systems. 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  4.1  describes  the 
family  of  plant  models  addressed  in  this  chapter  and  subsequently  details  some  key 
robust  stability  theorems.  The  main  results  o'  this  chapter  are  contained  in  Section  4.2, 
and  unstructured  multiplicative  uncertainty  versions  of  these  results  are  presented  in 


Section  :.i. 


4.1.  Preliminaries 


As  described  in  Chapter  2, plant  modeling  uncertainty  may  be  treated  by  representing 
the  true  physical  plant  by  a  family  Fg  of  possible  plant  models.  An  approach  (although 
possibly  conservative)  to  this  problem  is  to  allow  the  family  Fg  to  feature  dynamic 
(unstructured)  uncertainty.  Two  common  paradigms  for  dynamic  uncertainty  in  plant 
models  are  the  additive  and  multiplicative  representations  (see  Doyle  and  Stein  [  1 98 1  ]). 
For  instance,  additive  unstructured  uncertainty  is  represented  by  the  family  of  plant 
models  (2.2),  while  multiplicative  unstructured  uncertainty  is  represented  by  the  family 
of  plant  models 

Fg  =  |(I  +  A(z)  W,(z))  P0(z)  S  A(z)  II  <  5}  (4.1) 

along  with  identical  assumptions  as  in  (2.2). 

Kimura  [1984]  has  studied  the  problem  of  robustly  stabilizing  the  family  of  plants 
Fs  described  by  (2.2).  He  has  obtained  necessary  and  sufficient  conditions  for  robust 
stabilizability  of  a  family  of  SISO  plants  in  the  context  of  LT1  controllers.  (For  the 
MIMO  case,  see  Vidyasagar  and  Kimura  [1986].)  This  chapter  focuses  on  an  additive 
unstructured  family  of  plants  [somewhat  different  from  (2.2)]  of  the  form 

Fg  =  |P0(z)  +  AW(z):  II  A  II  <  8}  (4.2) 

where  P0(z)  and  W(z)  are  as  in  (2.2)  but  A  is  allowed  to  be  any  NLTV,  causal,  stable 
operator  with  II  A  II  <8. 

Remark  (4.1).  The  essential  difference  between  the  family  of  plant  models  (4.2)  and 
that  of  Kimura  [1984]  is  that  the  family  Fg  in  (4.2)  allows  for  nonlinear  perturbations  A, 
which  although  more  conservative,  may  be  a  reasonable  hypothesis  for  modeling  errors. 


(For  further  discussion  on  families  of  plants  featuring  nonlinear  perturbations  see 
Khargonekar  et  al.  [1987]  and  Poolla  and  ling  [1987].)  Also,  whereas  Kirnura  permits 
unstable  perturbations  (as  long  as  the  number  of  unstable  poles  remains  unaltered),  the 
perturbations  in  the  family  of  plant  models  (4.2)  are  restricted  to  be  stable.  Thus, 
strictly  speaking,  the  families  (2.2)  and  (4.2)  cannot  be  compared.  This  is  further 
discussed  in  Poolla  and  Ting  [1987].  It  shall  be  assumed  that  the  entire  family  of  plant 
models  (4.2)  is  strictly  proper.  This  assumption  is  merely  to  ensure  that  any  proper 
controller  one  designs  results  in  a  well-posed  feedback  system. 

□ 

The  first  step  in  treating  robust  stabilization  problems  of  families  of  plants  Fs  as  in 

(4.2)  is  to  obtain  robust  stabilization  theorems  (such  as  Theorem  (3.3))  for  families  of 
plants  featuring  nonlinear  perturbations.  The  following  result  is  a  straightforward 
generalization  of  the  results  of  Chen  and  Desoer  [1982]  to  the  family  of  plant  models 

(4.2) . 

Theorem  (4.1).  Consider  the  family  of  plants  F6  defined  by  (4.2).  Introduce  stable 
coprime  factorizations  of  P0  as  in  (3.7).  Then,  there  exists  a  (fixed)  LT1  controller  K  that 
robustly  stabilizes  F5  if  and  only  if 

inf  II W,XD  +  W2DQjDII  <  -  .  (4<3) 

Q,:LTI,  causa]  stable  5 

Proof.  The  necessity  follows  directly  from  Chen  and  Desoer  [1982]  by  noting  that  LT1 
perturbations  are.  in  particular,  N!.TV.  To  prove  the  sufficiency,  however,  a  little  care 
must  be  exercised  in  dealing  with  the  distributivity  of  nonlinear  operators:  suppose 
equation  (4.3)  holds.  It  then  follows  from  Adamjan  et  al.  [ 1 9781  that  there  exists  some 
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LTI.  causal,  stable  operator  Q,  which  achieves  the  infimum,  i.e.. 

II  W,XD  +  W,DQ,D  II  <  — .  Define  an  LTI  controller  K  =  (X  +  DQJIY-NQ,)"'  and  let 

"  ”  8  " 

T  =  K(1  +  P0K)-1.  From  Theorem  (3.1 )  it  is  clear  that  K  stabilizes  P0  and  using  (3.13b)  it 

follows  that  II  W,T  II  ^  .  From  Figure  3.1  it  is  easy  to  see  that  ej  =  u,— K(u,  +  Pe, ). 

8 

Since  K  and  (I  +  KP0)  1  are  linear,  this  becomes 

e,  =  (I  +  KP)-1  (u,-Ku,)  =  (I  +  TAW-,)-1  (1  +  KPor*  (u,-Ku2) 

=  W"1  (I  +  W,TA)_1  W,j(I  +  KP0)_1  u,  -  (1  +  KPor‘  Ku2}  . 

Notice  now  that  from  the  small-gain  theorem  (3.2),  (I  +  W2TA)  1  is  well  defined  and 
stable  for  all  A  with  II  A  II  <8.  Also,  since  W2,W2  \  (I  +  KP0)  1  and  (I  4-  KP0)  K  are  all 
stable,  one  may  conclude  that  the  operator  [d>,(P,K),  d>,(P,K)] :  [uj,u2]  -*e,  is  stable. 
Again,  from  Figure  3.1  since  P0  is  linear  it  follows  that 

e2  =  u2  +  Pe5  =  u2  +  P0ej  +  AW2e, 

=  u2  +  P0ut  —  P0Ke2  +  AW,e, 

=  (I  +  P0K)_1  (u2  +  P0u,  +  AW,ei) . 

T 

It  is  now  evident  that  since  K  stabilizes  P0,  the  operator  [<2>3(P,K),  d>4(P,K)]:  [u,,u2]  — *e2  is 
also  stable.  Thus,  d»(P.K)  is  stable  for  all  P  in  F5  (see  Definition  (3.2)),  and  this  concludes 
the  proof. 

□ 

Remark  (4.2).  For  NLTV  controllers,  the  proof  given  above  breaks  down.  The 
sufficiency  argument  relies  only  on  the  small-gain  theorem  and  remains  valid  even  for 
NLTV  controllers.  However,  the  necessity  argument  from  Chen  and  Desoer  [1982]  uses 


frequency  domain  techniques  and  is  no  longer  valid.  Indeed,  a  fundamental  result  of  this 
chapter  (see  Theorem  (4.4))  demonstrates  t ha t  the  small-gain  theorem  is  in  a  precise  sense 
both  necessary  and  sufficient  which  enables  one  to  prove  Theorem  (4.1)  for  NLTV 
controllers  as  well. 

□ 

4.2.  Main  Results 

In  this  section  it  is  shown  that  if  there  exists  a  NLTV  controller  that  robustly 
stabilizes  the  family  Fg  [see  (4.2)].  then  there  must  exist  a  linear  time-invariant 
controller  that  robustly  stabilizes  Fg.  A  similar  result  for  families  of  plants  featuring 
multiplicative  uncertainty  is  shown  in  the  next  section. 

For  purposes  of  clarity,  the  following  discussion  begins  by  outlining  the  sequence  of 
ideas  and  the  key  difficulty  in  the  proof  of  this  result.  To  this  end,  let  K  be  some  NLTV 
controller  that  robustly  stabilizes  the  family  (4.2).  Suppose  there  does  not  exist  any  LTI 
robustly  stabilizing  controller.  Then,  from  Theorem  (4. 1 ) 

inf  II  W,K(I  +  P0K)_I  II  >  -  . 

KiLTI.causal  8 

It  is  then  shown  (Theorem  (4.2))  that 

inf  II  W,K(I  +  P0Kr‘  II  =  inf  «  W2K(I  +  P0K)'' I  . 

K:NLTV.causal  ‘  K.LTI.causal 

-  -1  1 

Thus,  II  W,K(I  +  P0K)  II  > - be.  Now  by  applying  the  necessity  of  the  small-gain 

5 

theorem  (Theorem  (4.4))  it  is  possible  to  conclude  that  there  exists  a  causal,  stable 
perturbation  A  with  II  A  II  <5  and  such  that  (1  +  AW,K(]-rP0K )  1 )  1  is  unstable.  This  in 


turn  implies  that  K  does  not  robustly  stabilize  the  family  F5. 

Central  to  proving  this  main  result  is  the  necessity  of  the  small-gain  theorem 

(Theorem  (4.4)).  Given  a  (nontrivial)  NLTV  operator  M  with  II  M  II  >  —  +  e  it  is  quite 

5 

easy  to  construct  a  noncausal  operator  A  with  II  A  II  <8  and  such  that  (I+AM)  1  is 
unstable.  The  essential  difficulty  lies  in  exhibiting  a  causal  destabilizing  perturbation  A. 

In  order  to  prove  the  main  result,  Theorem  (4.5),  it  is  necessary  to  derive  several 
intermediate  results,  some  of  which  are  themselves  of  interest.  To  begin  one  may 
establish  the  following: 

Theorem  (4.2).  Let  A,  B,  and  C  be  causal,  stable,  linear  time-invariant  operators  with  C 
being  square  and  full  rank.  Let  8  =  618(5  and  C  =  CqC,  be  inner-outer  and  outer-inner 
factorizations  of  B  and  C,  respectively.  Then, 

inf  II  A  +  BQ,C  II  bh  =  inf  I  A  +  BQ2C  I  =  II  Bj*ACi*  II  H  .  (44) 

Q,  :  NLTV.causal.stable  2  Q2:LTI, causal, stable 

Proof.  The  proof  essentially  follows  from  arguments  due  to  Khargonekar  [1985].  First 

note  that  since  B,  and  (L  are  inner  (see  Arveson  [1975]),  Br  B,  =  I  and  C;  C|  =  I.  However, 

« 

except  in  the  trivial  case,  C.C;  ^  1  because  these  operators  are  defined  on  right-sided 
sequences.  Let  N  be  some  (fixed)  integer  and  define  the  operator  CN  =  C,  z  N.  Notice  that 
CN  is  a  contraction  and  define  the  operator  d»N  by 

C,CN  =  z_NI  +  d>N  .  (4.5) 

It  is  easy  to  see  that  lim  II  <f>N  II  =  0.  Let  Q,  be  any  NLTV,  causal,  stable  operator,  and 

N  — •  OO 

S 

notice  that  since  B  and  CN  are  contractions, 


II  A  +  BQjC  B  =  0  A  +  B,B0Q,C0Cj  II 

^  II  B,*ACN  +  B0Q,C0C,CN  II  . 

l  ei  Q2  :=  B^Cq  (which  is  a  causal,  stable  NLTV  operator).  Then,  from  (4.5), 

IIA  +  BQjCII  ^  IIB1‘aCn  +  Q,(z_NI +  (6n)II 

>  II  PnB,*ACn  +  PnQ2(z"NI  +  <2>n)  II 

=  I  PnB,*ACn  +  PNQ26N  II , 

where  the  last  equality  is  a  consequence  of  the  causality  of  Q2.  Continuing,  it  follows 


II  A  +  BQ,C  B  ^  II  P^'ac/z^  II  -  II  PNQ2d>N  II 

(4.6) 

>  #  PnB,*AC  ,VN  II  -  II  <2>n  II  K  , 

where  K  is  some  constant  independent  of  N.  It  is  easy  to  see  that 
II  PnB(  ACj  z  N  II  =  II  PnBi  AC,  (I  —  PN)  II .  Taking  limits  as  N  — ♦  oo  and  noticing  that 
lim  II  <&N  II  =  0,  one  immediately  obtains 

N"" *°° 

II  A  +  BQ,C  II  ^  lim  I  PnB,‘ac'(I  -  PN)  I  =  II  B,*AC,*  n  H  . 

N*“*oo 

*  * 

Since  Bj  .  A,  C,  .  and  PN  are  linear,  it  follows  that  there  exists  an  input  sequence  u  of 
arbitrary  norm  such  that  1!  ( A  +  BQ,C)u  U  /  II  u  II  ^  II  B,  AC,  II  H.  Thus, 


0,  :  NLTV. causal. s’.able 


"A  +  BQjCI*  >  H  B,  AC,  II  H 


In  point  of  fact,  it  is  shown  in  Adamjan  et  al.  [1978]  that 


A  +  BQ.C  II  =  II  B,  AC 


0,  :  LTl.causal.stable 


i  "  If 


Observing  that  LTI  operators  are,  in  particular,  NLTV  enables  one  to  conclude  equality  in 
(4.7),  proving  the  proposition. 


□ 

Remark  (4.3).  It  is  easy  to  sec  from  equation  (4.10)  that  given  any  e  >  0  and  any  real 
number  K  >  0,  there  exists  a  fixed  input  sequence  u  in  h2  such  that 


(A  +  BQC)u  II 
II  u  II 


B,  AC,  I  H  ~  e 


for  any  NLTV  causal  operator  Q  with  II  Q  II  <  K.  It  should  be  noted  that  a  result  similar 
to  Theorem  (4.2)  is  key  in  demonstrating  that  as  far  as  problems  of  uniformly  (or  H°°-) 
optimal  control  for  LTI  plants  are  concerned,  NLTV  controllers  offer  no  advantage  over 
LTI  controllers  (see  also  Khargonekar  and  Poolla  [1986]). 

□ 


The  next  step  is  to  show  that  operators  of  the  form  M  =  A+BQC  achieve  their  norm 
on  some  input  sequence  u  of  unbounded  energy.  This  is  demonstrated  in  the  following 
key  result. 


Theorem  (43).  Let  M  be  any  m -input,  p-output  NLTV  operator  of  the  form 
M  =  A  +  BQC  where  A,  B.  and  C  are  causal,  stable  linear  time-invariant  operators  with  C 
being  square,  full  rank  and  Q  is  any  causal,  incrementally  stable  NLTV  operator.  Let 
B  =  B,B0  and  C  =  C0C,  be  inner-outer  and  outer-inner  factorizations  of  B  and  C, 
respectively.  Suppose  further  that 


Bi’AC '  I!  „  >  |  +  €  • 
5 


(4.8) 


Then  there  exists  an  input  sequence  u  in  hJT  such  that 
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lim  II  PNu  II  =  oo  ,  and 

N— oo 


(4.9) 


lim 

N  — •  oo 


IIPNMull  , 


PNu 


>  -  +  €. 
8 


(4.10) 


Proof.  See  Appendix. 


□ 


Remark  (4.4)  .  Lssentially,  the  above  theorem  contains  a  "maximum  modulus  principle" 
for  a  certain  class  of  nonlinear  time-varying  operators,  i.e.,  it  states  that  the  norm  of  M 
is  achieved  on  the  "boundary"  of  h,  (i.e.,  h2e).  Recall  that  for  LTI  stable  operators,  the 
maximum  modulus  principle  asserts  that  the  norm  of  the  operator  is  achieved  on  some 
unbounded  sinusoidal  input  (i.e.,  the  axis).  In  point  of  fact  it  can  be  shown  that  the 
input  sequence  u  constructed  in  Theorem  (4.3)  can  be  guaranteed  to  be  periodic  .  Notice 

also  that  condition  (4.8)  implies  that  II  M  II  >  —  +  e  by  Theorem  (4.2). 

8 

□ 

The  above  result  is  used  to  prove  one  of  the  key  results  of  this  seciion.  In  particular, 
it  is  now  shown  that  the  small-gain  theorem  is  necessary  and  sufficient  for  nonlinear 
time-varying  operators  of  a  certain  form. 

Theorem  (4.4)  .  Let  M  be  any  NLTV  operator  with  II  M  II  >  —  as  in  Theorem  (4.3). 

8 

Then,  there  exists  a  NLTV,  strictly  causal,  stable  operator  A  with  II  A  II  <8,  and  such 
that  the  operator  (I  +  AM)-1  is  well-defined  and  unstable. 

Proof  .  It  follows  from  Theorem  (4.3)  that  there  exists  an  input  sequence  u  in  h such 
that  equations  (4.9)  and  (4.10)  hold.  This  in  turn  implies  that  there  exists  an  integer  N0 
such  that 


. 


II  PnMu  II 


(4.1  1) 


>  —  +  € .  N  >  Nn  . 

II  PNu  II  8 

Let  Mu  =  (y0,y,.y ,....)  =  y  and  define  the  sequence  u  =  -(0,0,0 . 0,uN(),uNo+, ,...).  For  a 

sequence  w  in  h2P,  define  the  integers  N(w)  =  smallest  integer  N  such  that  wN  ^  yN  ;  I(w) 
=  smallest  integer  i  such  that  w(  ^  0.  Define  the  operator  A  :  h,p  -*  h-,p  by 


A(w) 


PN<*->(U)  if  N(w)?±0 

z  A(z  w)  if  l(w)^0  . 


0 


otherwise 


(4.12) 


It  is  easy  to  verify  that  (4.12)  defines  a  strictly  proper,  stable,  nonlinear  time-invariant 

operator  A,  and,  it  follows  from  (4.12),  that  II  A  II  <  — 1 —  <  1.  Also,  (1  +  AM)"1  is 

8  +  e 

well-defined  (since  A  is  strictly  proper)  and 


(I  +  AM)u  =  u  +  Ay  =  u  +  u  =  (u0,u, . uN  ,,0,...) . 

Thus,  (I  +  AM)  (u0,u, . uNo_,,0,0,...)  =  u,  and  since  u  is  an  unbounded  sequence  (see 

equation  (4.9)),  it  follows  that  (I  +  AM)  is  an  unstable  operator.  This  completes  the 
proof. 

□ 

Remark  (4.5).  This  proposition  demonstrates  the  necessity  of  the  small-gain  theorem  for 
a  certain  class  of  NLTV  operators  (see  also  Vidyasagar  [1978]).  The  construction  of  the 
operator  A  above  basically  involves  determining  a  causal  inverse  of  M  on  the  input  u, 
which  exists  because  of  (4.11).  Essentially  this  construction  of  A  involves  "inverting"  M 
causally  on  this  input  u.  To  do  this  a  very  special  choice  of  the  input  u  must  be  selected 
in  order  to  ensure  that  II  A  II  <  8.  As  pointed  out  earlier,  determining  a  noncausal 


destabilizing  perturbation  A  is  quite  easy,  and  the  essential  problem  is  in  ensuring  that  A 


is  causal.  Ideally,  one  would  like  to  determine  a  linear  time-invariant  destabilizing 
perturbation  A.  This  problem  appears  at  first  glance  to  be  formidable  (see  also  Chapter  7). 


The  central  result  of  this  section  can  now  be  stated. 

Theorem  (4.5).  Consider  an  unstructured  additive  family  of  plants  F5  described  by 
(4.2).  Suppose  there  exists  a  NLTV  controller  K  that  incrementally  stabilizes  the  nominal 
plant  P0  and  robustly  stabilizes  Fs.  Then,  there  musl  also  exist  a  LT1  controller  K  that 
robustly  stabilizes  Fg. 

Proof.  (By  contradiction.)  Let  K  be  some  NLTV  controller  that  robustly  stabilizes  F6 
and  let  Q  be  its  corresponding  incrementally  stable  Youla  parameter  (see  Theorem  (3.1)). 
In  particular,  it  follows  that  K(I  -  Pl-K))-1  is  stable  for  all  P  in  F6  (see  Remark  (3.1)). 
Thus,  since  A  and  W  are  also  stable. 


xlf  =  I  +  AW(-K)(I-(P0  +  awx-k))' 


is  stable  for  all  A  with  II  A  II  <1.  Notice  now  that 


=  I  +  AW(— KXI  +  P0K)  1  (I  -  AW(— K )( 1  +  PoK)"1)' 


=  (I  —  AW(— K)(I  +  P0K)  ’)  1  . 


(4.13) 


Also,  from  (3.13b),  K(I  4- P0K)  —  XD  +  DQLX  where  D,  X,  and  D  are  obtained  from 
stable  coprime  factorizations  of  the  nominal  plant  P0.  With  A  = -WXD,  B  =  — WD. 
C  =  D  it  now  follows  from  (4.13)  that 


(1  —  AM)  is  stable  for  all  A  with  II  A  II  <  § 


(4.14) 


where  M  =  A  +  BQC . 


Now  suppose  that  there  does  not  exist  any  linear  time-invariant  controller  that 
robustly  stabilises  F.  This  implies,  from  Theorem  (3.3),  that 

1 

inf  II A  +  BQ,C  II  > - he.  (4.15) 

Q,  LTl  cdu«.al  stable  8 

However,  from  Adamjan  et  al.  [1978]  it  follows  that 

inf  II  A  +  BQ,C  II  =  II  b‘ac‘  II  H  . 

Q,  :  LTl.causal.stable 

«  »  1 

Therefore,  II  B,  AC,  II H  > - f  e,  and  consequently  from  Theorem  (4.4),  one  may 

8 

conclude  that  there  exists  a  causal,  stable  operator  A  with  II  A  II  <8  and  such  that 
(I  —  AM)  1  is  unstable.  This  renders  (4.14)  impossible,  concluding  the  proof. 

□ 

Remark  (4.6).  The  above  theorem  demonstrates  that  nonlinear  time-varying  controllers 
provide  no  advantage  over  linear  time-invariant  controllers  for  the  problem  of  robustly 
stabilizing  a  family  of  plants  featuring  dynamic  uncertainty  (see  (4.2)).  In  general, 
adaptive  controllers  are  nonlinear  and  time-varying;  therefore,  this  result  implies  that 
adaptive  controllers  are  not  useful  for  robustly  stabilizing  a  family  of  plants  featuring 
dynamic  uncertainty.  Due  to  the  proliferation  of  controller  synthesis  methods,  results  of 
this  type  are  invaluable  in  ruling  out  various  synthesis  methods  for  certain  classes  of 
control  problems. 

The  results  of  Theorem  (4.5)  are  not  entirely  unexpected  for  the  following  intuitive 
reason.  Due  to  the  nature  of  dynamic  uncertainty  (i.e.,  possibly  unbounded  McMillan 
degree  of  the  true  plant  model)  it  is  impossible  to  further  identify  the  true  plant  via 
input-output  experiments.  Hence,  a  nonlinear  time-varying  (or  adaptive)  control  strategy 
will  not  prove  advantageous. 


Remark  (4.7).  The  above  theorem  makes  critical  use  of  the  technical  restriction  that  the 
MLrV  robustly  stabilizing  controller  incrementally  stabilizes  the  nominal  plant.  This 
hypothesis  translates  to  stating  that  the  closed-loop  system  composed  of  the  NLTV 
controller  and  the  nominal  plant  be  "smooth,"  which  is  a  reasonable  requirement.  An 
interesting  open  question  to  examine  is  whether  or  not  Theorem  (4.5)  holds  for 
arbitrarily  NLTV  controllers. 

4.3.  The  Case  of  Multiplicative  Uncertainty 

This  section  extends  the  results  of  the  previous  section  to  families  of  plants  featuring 
unstructured  multiplicative  uncertainty.  Here,  the  family  of  plant  models  [somewhat 
different  from  (4.1  )]  is  of  the  form 

F8  =  (( 1  +  A)P0(z):  II  All  <  8}  (4.16) 

where  P0(z)  and  A  are  as  in  (4.2). 

For  the  family  of  planls  ( 4. 1 )  featuring  multiplicative  unstructured  uncertainty,  one 
may  derive  the  following  robust  stability  criteria  analogous  to  the  results  of  Theorem 
(4.1). 

Theorem  (4.6).  Consider  the  family  of  plants  F5  defined  by  (4.1).  Introduce  stable 
coprime  factorizations  of  P0  as  in  (3.7).  Then  there  exists  a  (fixed)  LTI  controller  K  that 
robustly  stabilizes  F5  if  and  only  if 

-  -  1 

inf  II  NX  +  NQ,D  II  <  -  . 
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The  proof  of  'his  result  is  very  similar  to  that  of  Theorem  (4.1 )  and  is  hence  omitted. 
Again,  however,  the  pnxtf  breaks  down  in  the  case  of  NLTV'  controllers  (sec  Remark 


(4.2)),  and  it  is  necessary  to  employ  the  results  of  Theorem  (4.4)  in  order  to  generalize 
ihese  results  to  NLTV  controllers.  The  following  necessary  condition  for  the  robust 
stabilizability  of  a  family  of  plants  Fg  as  described  by  (4.16)  can  now  be  stated. 

Theorem  (4.7).  Consider  the  family  of  plants  Fg  described  by  (4.16)  and  let  the  nominal 
plant  P()  admit  a  stable  coprime  factorization  as  in  (3.7).  Let 
M  =  — P0K(I+P0K)  1  =  — N(X+QD)  be  expressed  linearly  in  Q  as  M  =  A+BQC  by  setting 
A  =  —NX,  B  =  — N.  and  C  =  D.  Then  there  exists  a  NLTV  controller  K  that  robustly 
stabilizes  Fg  only  if  the  operator  C(l— AM)  1  is  stable  for  all  A  with  II  A  II  <5. 

Proof.  Suppose  K  is  a  NLTV  operator  that  robustly  stabilizes  the  family  of  plants  Fg. 
From  Remark  (3.1)  it  follows  that  for  any  plant  P  in  Fs  the  operators  ( I — P( — K ))  1  and 
K(I— P(— K))  '  are  both  stable.  It  is  straightforward  to  express 

(I  -  P(-K))-1  =  (I  -  (I  +  A)  P0(— K)f 1  =  (l  +  P0K  -  AP0(-K))"' 

=  (l  +  P0Kr’  ( I  —  AM)-1 

K(I  -  P(-K))“‘  =  K(1  +  P0K )-1  (I  -  AM)"'  . 

Employing  (3.13a)  and  (3.13b)  and  noting  that  the  addition  and  composition  of 
stable  operators  is  stable  it  follows  that 

NKd-P(-K))"1  +  D(I— P(— K))_I  =  [N(X  +  DQ)  6  +  D(Y-NQ)  6]  (I-AM)-1 

=  IX  I-AM)-1 

is  indeed  a  stable  operator  for  all  A  with  II  A  II  <5.  This  completes  the  proof. 

□ 

One  may  now  derive  the  following  result. 


Theorem  (4.8).  Consider  the  problem  of  robustly  stabilizing  the  family  of  plants  Fg 


described  in  (4.16).  Suppose  there  exists  a  NLTV  controller  K  that  incrementally 

stabilizes  the  nominal  plant  P0  and  robustly  stabilizes  Fg.  Then,  there  must  also  exist  a 

LTI  controller  K  that  robustly  stabilizes  Fs. 

Proof.  (By  Contradiction).  Let  the  operator  M  =  A+BQC  be  as  defined  in  Theorem 
(4.7).  Let  K  be  some  NLTV  controller  that  robustly  stabilizes  Fg  and  let  Q  be  its 
corresponding  incrementally  stable  Youla  parameter  (see  Theorem  (3.1)).  From 
Theorem  (4.6)  this  implies  that  the  operator  C(I— AM)  1  is  stable  for  all  A  with  II  A  ii  <8. 

Now  assume  that  there  does  not  exist  any  LTI  controller  K  that  robustly  stabilizes 
Fg.  This  implies,  from  Theorem  (3.3),  that  for  some  fixed  c>0 

*  «  1 

inf  II  A  +  BQC  II  =  B  Bj  AC,  B  H  >  —  +  €  . 
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Consequently,  from  Theorem  (4.3)  one  may  conclude  that  there  exists  an  input  sequence 
u  in  h;™  satisfying  (4.9)  and  (4.10)  and  such  that  Mu  =  y.  Following  Theorem  (4.4)  it  is 

possible  to  define  a  stable  causal  operator  A  similar  to  (4.12)  such  that  II  A  B  <  —  <  — 

8+c  8 

and  (I— AM)  1  is  well-dehned  and  unstable. 

The  operator  (I— AM)  1  maps  (I— AM)  '(u0,u, . uN  _,,0,0....)  =  u  and  thus 

C(  1 — AM)  1  (u0,u, . uN  ,,0.0,...)  =  Cu.  Now  suppose  that  the  sequence  Cu  is  bounded. 

It  is  known  that 

Mu  =  (A  +  BQC)u  =  Au  +  BQCu  =  -NXu  +  BQCu 
=  — (1— YD)u  +  BQCu  =  -u  +  YDu  +  BQDu  . 

Thus,  if  Cu  =  Du  is  bounded  then  YDu  +  BQDu  is  bounded  since  Y.  B,  and  Q  are  ail  stable 
operators.  Hence, 
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Mu  =  — u  +  v  where  v  €  h,  . 


ijj 

IS 


I  his  implies  that  II  Mu  II  /  II  u  II  =1  which  contradicts  the  fact  that  the  input  sequence  u 

II  Mu  II  1 

was  chosen  such  that  -  > - he.  Therefore,  the  sequence  Cu  cannot  be  bounded. 

II  u  II  8 

Now  assume  that  the  sequence  Cu  is  unbounded.  Then  the  operator  C(I— AM)  1 
which  maps  C(I— AM)  1  (u0,u, . uN  ,0,0,...)  =  Cu  is  an  unstable  operator.  From 


Theorem  (4.6)  it  follows  that  there  does  not  exist  any  NLTV  robustly  stabilizing 
controller  for  F5  and  this  concludes  the  proof. 


CHAPTER  5 


ROBUST  STABILIZATION  PROBLEMS 
WITH  ADDITIVE  BLOCK-STRUCTURED  UNCERTAINTY 

The  research  in  this  chapter  is  focused  upon  various  aspects  of  the  robust 
stabilization  problem  with  additive  block-structured  uncertainty.  As  is  well-known,  the 
solution  of  this  problem  reduces  to  an  optimization  problem  of  the  form 

inf  il  A(A +  BQC)  A-1  I!  „  .  ,  . 

A.A'1,  Q  in  Hoj 

Here,  A  is  constrained  to  be  block-diagonal  conformally  to  the  uncertainty  structure  and 
A,  B,  and  C  are  determined  from  a  stable  coprime  factorization  of  the  nominal  plant 
model.  An  approach  to  the  solution  of  this  infinite  dimensional  optimization  problem  is 
the  so-called  A-Q  iteration  of  Doyle  [1983].  The  basic  idea  involved  here  is  to  alternately 
freeze  the  scaling  A,  find  the  best  Q,  freeze  Q,  find  the  best  A  and  so  on.  While  the  above 
problem  is  convex  individually  in  A  and  in  Q, it  fails  to  be  jointly  convex.  Hence  there  is 
no  guarantee  that  this  procedure  will  converge  to  a  global  minimum.  Further,  A-Q 
iteration  has  some  numerical  drawbacks.  For  instance,  at  each  step,  the  optimal  scaling 
A(jg>)  must  be  computed  at  each  frequency.  Also,  A(jo>)  must  then  be  approximated 
rationally  in  order  to  determine  the  optimum  Q. 

In  this  chapter  the  infinite  dimensional  optimization  problem  presented  above  is 
reduced  to  a  nonlinear  optimization  problem  involving  a  finite  number  of  variables.  This 
reduction  is  accomplished  by  exploiting  a  key  observation  due  to  Safonov  [1986]  that  the 
optimization  problem  in  (5.1)  is  dependent  only  on  the  values  of  A  and  A-1  at  the 
unstable  poles  of  the  nominal  plant  model  P0.  Using  this  fact  it  is  possible  to  significantly 
reduce  the  computational  effort  involved  in  this  controller  design  problem.  Also,  this 


approach  provides  insight  into  the  robust  stabilization  problem  with  additive  block- 
structured  uncertainty.  In  particular,  some  simple  sufficient  conditions  on  the  nominal 
plant  P0  are  derived  which  guarantee  infinite  additive  stability  margin. 

In  a  receni  paper,  Safonov  [1986]  has  treated  a  more  general  optimization  problem  of 
which  robust  stabilization  with  block-structured  uncertainty  is  a  special  case.  He  has 
been  able  to  show  that  this  more  general  problem  also  reduces  to  nonlinear  programming 
with  a  finite  number  of  variables.  However,  when  these  results  are  specialized  to  the 
robust  stabilization  problem  treated  in  this  chapter,  the  explicit  dependence  of  this 
optimization  problem  on  the  scaling  function  residues  is  complex.  Thus,  one  of  the 
contributions  of  this  chapter  may  be  viewed  as  a  simplification  of  the  particular  nonlinear 
programming  problem  that  must  be  solved  for  robust  stabilization  with  block-structured 
uncertainty. 

A  key  step  in  recovering  the  controller  designed  via  the  techniques  of  this  chapter  is 
the  determination  of  a  stable,  minimum  phase  (outer)  function  subject  to  certain 
interpolation  constraints  at  the  unstable  poles  of  P0.  The  degree  of  the  resulting 
controller  is  critically  dependent  on  the  order  of  these  outer  functions.  A  new  method  for 
computing  these  interpolants  that  yields  far  lower-order  interpolants  than  the  method  of 
Youla  et  al.  [1974]  shall  be  presented. 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  5.1  contains  a  detailed 
account  of  the  residue  iteration  algorithm  and  some  remarks  on  the  numerical  aspects  of 
this  problem.  The  case  of  infinite  stability  margin  is  examined  in  Section  5.2.  Following 
this,  in  Section  5.3  these  results  are  illustrated  via  a  few  numerical  examples. 


5.1.  Main  Results 


It  will  be  shown  in  this  section  that  the  problem  of  robustly  stabilizing  an  additive 
structured  family  of  square  plants  reduces  to  a  computationally  tractable  linear  algebraic 
optimization  problem  referred  to  as  residue  iteration.  As  pointed  out  in  Chapter  3,  many 
families  of  plant  models  (possibly  non-square)  with  block-structured  uncertainties  can  be 
transformed  into  equivalent  fictitious  square  families  of  plant  models  with  block- 
structured  uncertainties.  Some  numerical  aspects  of  this  problem  shall  also  be  addressed. 
As  stated  earlier,  attention  shall  be  focused  on  the  case  of  three  or  fewer  blocks  in  order 
that  the  results  obtained  are  both  necessary  and  sufficient.  Also,  the  following  exposition 

shall  deal  only  with  scalar  block  sizes  (i.e..  r—l  for  i-1,2 . k)  to  avoid  cumbersome 

notation  though  there  is  no  essential  loss  of  generality  by  doing  so.  For  additional 
information  on  the  results  presented  in  this  chapter  see  Ting.  Cusumano.  and  Poolla 
[1987], 


5.1.1.  The  Residue  Iteration  Algorithm 

Beginning  with  the  alternate  characterization  of  robust  stabilizability  as  in  Theorem 
(3.4)  the  following  result  is  easily  obtained. 

Theorem  (5.1).  Consider  the  family  of  plants  (2.3)  with  block-structured  uncertainty. 
Then,  there  exists  a  controller  K(z)  that  robustly  stabilizes  the  family  (2.3)  if  and  only  if 

sup  g{ r(AP0W  1  A  1 ))  >  8  .  (5  '’l 

a  e  Q 

Proof:  From  Theorem  (3.5)  it  follows  that  there  exists  a  robustly  stabilizing  controller 
if  and  only  if  inequality  (3.18)  holds.  Define  the  transfer  function  PN(s)  =  AP0A_I.  It 
can  be  verified  that  a  stably  coprime  factorization  of  PA  is 


PA  =  ( ANA  1 KADA  ')  1  =  ( ADA  ')  ‘(ANA  ') 
(ANA-‘)(AXA_I)  +  (ADA-1  )(AYA_1)  =  1  . 


It  is  now  evident  from  Theorem  (3.3)  that  inequality  (3.18)  holds  if  and  only  if  the 
unstructured  family  of  plant  models  |PX  4-  AW  :  II  A  II  ^  8}  is  robustly  stabilizable.  Now 
from  Theorem  (3.4)  this  family  is  robustly  stabilizable  if  and  only  if  inequality  (5.2) 
holds,  completing  the  proof. 

□ 


It  is  crucial  to  note  that  T(AP0W  *A  ’)  depends  only  on  the  residues  of  A  at  the 
unstable  poles  of  P0.  The  following  key  proposition  demonstrates  the  relationship 
between  the  smallest  Hankel  singular  value  of  AP0\V  '.\  1  and  the  residues  of  A  at  the 
unstable  poles  of  P0. 

Proposition  (5.1).  Let  ~(A.B,CO)  be  a  canonical  realization  of  (P0W  1 )+  with  A  in 
Jordan  form  and  let  the  McMillan  degree  of  (P0W  ')+  =  N.  Let 

A(s)  =  diag(dj(s),d,(s) . d^s))  and  for  i=l,2,...,k  define  the  NxN  residue  matrices  R; 

by 


Rj  =  diag{dj(p1),di(p2) . dt(  pN ) } 

where  the  diagonal  elements  of  A  are  p,  ,p2 . pN.  Tor  i-1,2 

B  and  let  c;  =  jth  row  of  C.  and  solve  the  Lyapunov  equations 

AM;  +  M,AT  =  bp? 

ATW;  +  W;A  =  ciTc, 


k  let  b;  =  ith  column  of 


(5.3) 


f or  the  grammians  M(  and  W(.  Define 
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M(R)  =  X  Rj  '  M,R,  ’ 


W(R)=  £  R  \VjR|  . 


(5.4) 


Then. 


o[r(AP0W  'A  '))  =  A.'*[M(R)W(R)]  . 


Proof.  Let  e,  =  i  column  of  the  kxk  identity  matrix  lh.  With  the  above  (tedious*) 
notation,  it  follows  from  elementary  matrix  algebra  that  a  realization  for  (AP0W-IA-1)+ 
is  rr(  A,BvCa.O)  where 


BA  =  Z  Rf '  bftT  CA  =  £  e,c,R, 


Also  notice  that  since  { e,  :  i  =  1,2 . k}  is  an  orthonormal  set. 


baBa  -  Z  Rf '  b,b,T  R  .  CatCa  =  £  R,c,t  c,R, 


Let  M(R)  and  W(R)  be  solutions  of  the  Lyapunov  equations 


(5.5) 


AM(R)  +  M(R)AT  =  BaBa 
ATW(R)  +  W(R)A  =  CaT  Ca  . 


(5.6) 


Then  from  well-known  results  (see  Glover  [1986]).  o.(r(AP0W  A  ))  =  A.  (M(R)W(R)). 
It  remains  to  be  shown  that  the  solutions  of  the  Lyapunov  equations  (5.6)  are  given  by 
the  equations  (5.4).  The  demonstration  of  this  fact  rests  on  the  key  observation  that  for 

i  =  1,2 . k,  A  and  R,  commute.  To  see  this,  let  J,  ,  J2 . J,  be  the  Jordan  blocks  of  A 

and  partition  R,  conformally  with  these  blocks  as 


Notice  now  that  R1  —  di(p])Ir  where  is  the  pole  corresponding  to  the  r-Xr^  Jordan  block 
Jj.  1  hus,  clearly  ARj  =  R;A.  Now  to  complete  the  proof  observe  that 


AT(£  R,WiRj)  +  (£  RWjRj)A 


=  I  R,(ATW,  +  WjA)R, 


=  I  =  Cat  Ca  . 


Thus,  the  solution  W(R)  to  the  Lyapunov  equation  (5.6)  is  given  by  W(R)  =  £RjWiRj- 


One  can  similarly  derive  the  second  part  of  equation  (5.4)  completing  the  proof. 


The  above  proposition  suggests  an  algorithm  for  solving  the  robust  stabilization 
problem  with  additive  block-structured  uncertainty  which  is  outlined  below. 

Consider  the  family  of  plant  models  (2.3). 

Step  I.  Compute  a  canonical  realization  tt(A,B,C,0)  for  (P0W-1)+  with  A  in  Jordan 


f orm.  say , 


Step  II.  For  i— 1,2 . k.  let  =  ith  column  of  B  and  let  c,  =  ith  row  of  C  and  solve  the 

Lyapunov  equations  (5.3)  for  the  grammians  M,  and  Wr 
Step  III.  Compute 


sup  \  (M(R)W(R))  =  y 
V '  • 


(5.7) 


where  M(R)  and  W(R)  are  as  in  (5.4).  Let  Rj° . RK°  be  a  solution  to  this  optimization 

problem. 

Step  IV.  For  i=L2 . k  find  outer  (stable,  minimum  phase)  functions  d;(s)  that 

interpolate  p,  to  R,°(j,j).  Such  interpolants  always  exist  (see  Youla  et  al.  [1974]).  Define 
A(s)  =  diag|d,(s) . dk(s)}. 

Step  V.  Let  PU  =  (AP0W  *A  ')+  and  let  Ps  =  AP0W  ]A  Pu.  Determine  via  the 
methods  of  Glover  [1984]  an  optimal  N-lst  order  anticausal  Hankel  norm  approximant 
PN_,  for  Pu  such  that  Pu— PN_,  is  inner,  where  N  -  McMillan  degree  of  Pu. 

Step  VI.  Then  an  optimal  robustly  stabilizing  controller  for  the  family  (2.3)  is  given  by 


K(s)  =  -A"'(Pn_i  +  Ps)  'a 


(5.8) 


The  above  residue  iteration  algorithm  is  illustrated  in  Section  5.4  via  a  few  examples. 
The  development  now  shifts  to  address  some  important  practical  details  concerning  this 
algorithm. 


5.1.2.  Numerical  Aspects  of  Residue  Iteration 


The  key  numerical  bottleneck  in  the  above  algorithm  is.  of  course,  Step  III.  In 
general  it  is  not  true  that  the  functional  A_  (M(R)W(R))  is  concave  in  the  parameter 

vector  R  =  (Rj.R, . Rk).  Trivial  counterexamples  to  demonstrate  this  can  be  easily 

generated.  This  algorithm  has  been  implemented  using  a  multidimensional  line  search 
without  derivatives  such  as  the  method  of  Hooke  and  Jeeves  (see  Bazaraa  and  Shetty 
[1979]).  Several  numerical  examples  indicate  that  the  functional  A./4(M(R)W(R))  may. 
however,  be  quasiconcave  in  R  (see  Rockafellar  [1970]).  This  would  be  very  nice  since 
quasiconcave  functions  have  at  most  one  local  maximum.  Thus,  upon  arrival  at  a  local 
maximum  in  Step  III  one  need  only  check  the  boundary  behavior  of  \/‘(M(R)W(R))  as  R 
approaches  00  to  determine  the  global  maximum.  These  issues  are  still  under 
investigation. 

In  general, the  residues  d^)  i  =  1,2 . k  ,  j  =  1,2 . N  will  be  complex.  For 

real  poles  p-  it  suffices  to  carry  out  the  maximization  (5.7)  over  real  positive  residues.  It 
should  also  be  noted  that  as  in  Doyle  [1982]  one  may,  without  loss  of  generality,  choose 
d,(s)  =  1. 

Step  IV  of  the  above  algorithm  is  the  next  key  step  and  this  involves  determining 
outer  interpolating  functions.  The  degree  of  these  functions  directly  affects  the  McMillan 
degree. of  the  controller  (in  Step  V!)  and,  therefore, it  is  desirable  to  have  low-order 
interpolants.  The  interpolation  technique  of  Youla  et  al.  [1974]  is  unsatisfactory  in  this 
regard.  Also,  the  C  or  logarithmic  interpolation  methods  of  Bali  and  Helton  [1979], 
which  provide  minimum  norm  outer  interpolants,  yields  transcendental  functions  which 
then  have  tc  be  rationally  approximated.  The  following  exposition  outlines  a  technique  to 
determine  outer  interpolating  functions  that  always  yield  (much)  lower-order  rational 


outer  functions  than  the  methods  of  Youla  and  is  computationally  much  less  demanding 
than  the  methods  of  Ball  and  Helton. 

Let  [a,  :  i  =  1.2 . L}  and  [b, :  i  =  1.2 . L)  be  two  conjugate  closed  sets  of  complex 

numbers  in  H  and  C-{0[,  respectively.  Suppose  further  that  if  a,  is  real,  then  bj  is  in  R+ 
for  the  corresponding  integer  i.  Then  it  is  well-known  (see  Youla  et  al.  [1974])  that  there 
exists  an  outer  (stable,  minimum  phase)  function  f(s)  such  that 

f(ai)  =  bi  i  =  1.2 . L.  (5.9) 

For  any  positive  integer  in  m  define  the  Nevanlinna-Pick  matrix 

,  l/m  ,  r  I/m 
bj  4-  b- 

- -  i.J  =  1 .2 . L 

a,  +  aJ 


where  a  denotes  complex  conjugate  of  a  and  bWm  is  the  principal  m,h  root  of  b.  It  is  now 
possible  to  show  the  following  result. 

Theorem  (5.2).  Consider  the  problem  of  outer  function  interpolation  described  above  in 
(5.9).  Let  m0  be  the  smallest  integer  such  that  N  is  positive  definite.  Then  there  exists 

o 

an  outer  rational  function  f(s)  of  order  at  most  m0L  that  solves  this  interpolation 
problem.  One  such  interpolant  is  given  by 


Here  g(z)  is  determined  by 
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L  fix 

t  'l  l 

(— 1 )  ( z - ) 


g(z)  = 


.-i  z_a. 


7T  (  I  Qfj  i  -  £ 


i=.  Z“a. 


(5.10) 


1  -  fi.fi, 


1  —  or  .Of. 

1  j 


(  7T  I  I  ) 
1=1 


1/m 

ai —  1  Q  °-l 

<*,  = -  .  fij  = -  .  i  =  1 ,2 . L 

a‘+1  b 


Proof:  First  note  that  since 


lim  N_  = 

m 


>  o, 


!  a,  +  a. 


there  must  exist  a  smallest  integer  m0  such  that  Nm  >0.  Defining  and  fi|t  i  =  1,2 . L 

as  in  (5.10),  it  is  easy  to  see  using  a  H— *D  conformal  map  that  N  >0  implies  that 


1  -  fi.fij 

M  =  - -  >  0  . 


(5.1 1) 


1  —  or  Of 

i  j 


This  in  turn  implies  via  Krein  and  Nudelman  [1977]  that  for  some  e>0  there  exists  an 


analytic  function  g(z)  :  D  — ♦  D,_e.  Mere  D,_e  is  the  open  disk  of  radius  1  — c  in  C.  This 
result  makes  critical  use  of  the  positive  definiteness  of  the  matrix  M  (see  (5.11)).  One 
such  rational  function  g(z)  of  order  at  most  L  may  he  determined  using  formulae  (5.10). 


Notice  now  that 


,  s-1 

1  +  g( - ) 

.  ,A  s+I 

h(s)  =s -  :  U  — *  IT 

1  -g(  — ) 
s+1 


is  an  outer  function  interpolating  h(ai)  =  bI  °  for  i  =  1,2 . L.  It  is  now  evident  that 

mo 

f(s)  =  h(s)  solves  the  outer  function  interpolation  problem  (5.9)  and  is  a  rational 
function  of  order  at  most  m0L. 

Example  (5.1).  Consider  the  problem  of  finding  an  outer  f unction  f(s)  such  that  f(  1 )  =  1. 
f(2)  =  3,  f  (3)  =  2.  It  can  be  verified  that  the  method  of  Youla  el  al.  [1974]  yields  a  48th 
order  rational  function.  However,  using  the  method  of  Theorem  (5.2)  one  obtains 


N_  = 

m 


3  m  +2  m 


3  m  +2  m 


and  N'j^O.  Thus,  it  is  possible  to  compute  a  15th-order  outer  rational  function  f(s)  that 
meets  the  required  interpolation  constraints. 


5.2.  The  Case  of  Infinite  Stability  Margin 


In  this  subsection  some  sufficient  conditions  on  the  nominal  plant  P0  will  be  derived 
in  order  that  it  exhibit  infinite  additive  block-structured  stability  margin.  These  results 
hold  for  an  arbitrary  number  of  blocks  (even  for  ^  3  blocks).  In  particular,  one  may 
demonstrate  the  following  interesting  result: 

Theorem  (5.3).  Let  P0  be  a  transfer-function  matrix  with  McMillan  degree  (P0)+ =  N. 

and  consider  an  additive  block-structured  family  of  plant  models  as  in  (2.3).  Partition  P0 

conformally  with  its  block  diagonal  uncertainty  structure,  let  D0  =  the  block  diagonal 

A 

part  of  P0,  and  define  Q0  =P0— D0.  Suppose  further  that  the  following  technical  condition 
holds: 


For  i^i,  ,  and  j^j,  the  unstable  poles  of  Q0(i,j)  and  the  unstable 
poles  of  Q0(i,,j, )  are  distinct.  If  McMillan  degree  of  (Q0)+  =  N  (5.12) 

then  the  maximal  achievable  stability  margin  is  infinite. 

Proof.  Without  loss  of  generality  assume  W(s)  =  1  in  (2.3).  As  a  consequence  of  the 
technical  hypothesis  (5.12),  it  follows  that  there  exist  outer  functions  A  in  fi  with 
arbitrary  residues  d2(pj)  ==  R2( j, j )•  In  particular,  for  any  fixed  real  number  a  there  exists 
an  outer  function  An  in  Q  such  that 

(w;,).=(D„+«QA. 

It  now  follows  from  Glover  [1984]  that 

yo  =  <L(r(A„p0A;1))=  inf  nr(AQp0A;1-pN_1)ii  . 

Pn-i  =  of 
McMilUn  dffiff  N-l 

Suppose  this  infemum  is  achieved  by  P^_,.  Then, 
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y0=  «nAap0Aft1-pN°_1)ii  =  iir(D0  +  aOo-PN0-.) 


^  II  r(aQ0-P°_,)ll  -  II  r(D0) 


inf  II  r(aO„-PN_,)ll  -  II  r(D0)ll 

PN-I  :  °f 
McMillan  deftet'N-1 


=  a  o(Q0)-  II  r(D0) 


where  the  last  step  above  critically  uses  the  hypothesis  that  the  McMillan  degree 
(Q0)+  =  N.  Since  o.(r(Q0))  ^  0  it  follows  that 


y  =  sup  o.(  AP0A  ' )  ^  sup  yn  =  oo 

A  €  fi  n 


or  that  (by  Theorem  (5.1))  the  maximal  achievable  stability  margin  is  infinite. 


□ 


Essentially  what  the  above  theorem  states  is  lhai  if  every  "unstable  mode"  of  P0  is 
present  in  the  off-diagonal  blocks  of  P0  then  the  maximal  achievable  stability  margin  is 
infinite.  It  is  possible  to  prove  much  sharper  forms  of  Theorem  (5.3)  (by  weakening  the 
technical  hypothesis  (5.12));  however,  these  results  involve  cumbersome  notation  and  are 
not  as  clean  as  the  above  theorem  (see  Cusumano  [1987]). 

Example  (5.2).  Consider  the  family  of  plant  models  as  in  (2.3)  with 


po  = 


(2s— 5) 


15 


( s — 2 )( s — 3 )  s— 3 


15 
s— 2 


1 

s+ 1 
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Here.  Q  =  15 


and  McMillan  degree  (P0)+  =  McMillan  degree  (Q0)+  =  2  and 


from  Theorem  (5.3)  it  is  easy  to  see  that  P0  can  be  robustly  stabilized  against  arbitrarily 
large  additive  block-structured  perturbations. 


5.3.  Numerical  Examples 

This  section  contains  a  few  numerical  examples  to  illustrate  some  of  the  concepts  of 


the  previous  sections. 

Example  (5.3).  Consider  the  2-input  2-output  transfer  function  matrix 


s—  1  s— 2 


s— 3  s— 2 


Consider  also  the  additive  family  of  plant  models  with  structured  uncertainty  (as  in 


(2.3)): 


Ai  0 

';s=  l’o+  0  4t  : 


Suppose  the  objective  is  to  find  the  largest  stability  margin  5  against  which  robust 
stabilization  is  possible  using  the  residue  iteration  algorithm  of  Section  5.2. 

Step  I:  A  canonical  realization  for  (P0)+  with  A  in  Jordan  iorm  is  jt(A,B,C.D)  where 


•  *  .  j  ’j  '  *»  A-  ’  •  •  • 


'0  0  i 


U  ni 


10  .5  1 


10  0) 


Step  II:  Solving  the  Lyapunov  equations  (5.3)  for  the  grammians  M1  and  W.  one  obtains 


.5 

0 

.5 

0 

0 

0 

0 

0 

0 

m2=  (o 

.25 

0 

.5 

0 

.67 

0 

0 

0 

.5 

.33 

0 

0 

0 

0 

.33 

.25 

0 

w2=  0 

1 

.4 

0 

0 

0 

0 

.4 

.17 

Step  III:  Solving  the  optimization  problem  (5.7)  yields 

*  0  0 

R.  =  I.  and  R,  =  0  .952  0  . 

b  0  .85 

Also,  the  supremum  y  =  0.123. 

Step  IV:  Using  the  method  outlined  in  Section  5.2  it  is  easy  to  see  that  the  Nevannlina- 
Pick  matrix  Nm  associated  with  this  interpolation  problem  is 


(.952) m  -f  (.85) ' 
5 


(.952) m  (.952) m  +  (.85) m 

N  =  2  ^ 

Nm  1  1  1 

(.952)™  -f  (.85)  m  ( .85) m 

5  3 

Further,  N,  is  positive  definite  and  using  the  methods  of  Krein  and  Nudelman  [1979]  it  is 
possible  to  find  a  second-order  rational  outer  interpolant.  One  such  f unction  is 

.5 1 03s2  +  . 5746s +  4.065 


,8457s  +  .3614s  +  L515 


^Swl'i'wivW1 
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A(s)  = 


0 

d(s) 


Step  V:  Using  the  methods  of  Glover  [1984]  an  optimal  second-order  anti-causal  Hankel 
norm  approximant  P,(s)  of  ( AP0A-1 )+  may  lx;  found  to  be 


P,(s)  = 


.1 149(s*  +  21.85s  +  552.7)  -,043(s~  -  436s  +  4270) 


.-.043(s2  -  257.9s  +  1042)  -.1 149Cs2  -  24.7s  +  658.8) 


(s  —  80.4)  (s  —  1.31) 


Step  VI:  A  realization  for  an  optimal  robustly  stabilizing  controller  can  be  found  via 
(5.8)  to  be  K(s)  =  ;7(A,B,C,D)  where 
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Example  (5.4).  With  P0  in  Example  (5.3  )  consider  the  structured  family  of  plant  models 


H8,.8,)  = 


!i>  *+* 

A, 

G 

o 

0 

^2 

<8: 


This  example  shall  demonstrate  the  tradeoffs  between  the  achievable  stability  margins  8, 
and  8,  in  the  two  channels.  This  can  be  done  by  a  rescaling  of  the  nominal  plant  by  the 

ratio  q  =  —  and  then  using  the  methods  of  Section  5.3.  Since  the  details  involved  are 

S, 

identical  to  those  in  the  previous  example  only  the  final  results  shall  be  discussed. 
Plotted  in  Figure  5.1  is  a  curve  a  (shown  in  the  solid  line)  which  represents  the  maximum 
achievable  stability  margins  on  each  channel  as  a  function  of  q. 


Figure  5. 1 . 


Maximum  achievable  stability  margins  on  each  channel. 


Any  desired  stability  margin  below  this  curve  can  be  obtained  by  some  FDLTI  controller 
K(s).  Following  the  works  of  Khargonekar  el  al.  [1987]  and  Pool  la  and  Ting  [1987],  an 
interesting  open  issue  is  the  possible  advantage  of  the  use  of  nonlinear  time-varying 
control  towards  obtaining  better  stability  margins.  Suppose  one  finds  a  perturbation  pair 
(A^.A,0)  with  norms  (  II  A,°  II  ,  H  A,°  II  )  =  (81°,8,0)  such  that  P0  +  diagl  Aj.A-,}  has  one  less 
unstable  pole  than  P0.  Then,  it  follows  from  an  elementary  application  of  the  arguments 
of  Khargonekar  et  al.  [1987]  that  there  cannot  exist  any  (possibly)  nonlinear  time- 
varying  controller  that  yields  stability  margins  of  (8, ,8,)  on  the  two  channels  where 
8,  ^  8,°  and  8,  ^  8,°.  Using  ad  hoc  techniques  it  can  be  shown  that  for  P0  as  above  there 
exist  perturbation  pairs  (A,.A2)  of  norms  (0.25,  0);  (0.167,  0.05)  and  (0.08,  0.25)  such 
that  P0  +  diag( A, .A,)  has  one  less  unstable  pole  than  P0.  Also,  an  (exhaustive!)  search 
reveals  that  these  are  the  "smallest"  size  diagonal  perturbations  that  decrease  the  number 
of  unstable  poles  of  P0  (in  terms  of  the  partial  ordering  on  R^).  Thus,  one  may  conclude 
that  for  desired  stability  margins  above  the  curve  /3  (shown  in  the  dotted  line)  there  do 
not  exist  any  controllers  nonlinear  and/or  time  varying  that  achieve  these  margins. 
However,  for  desired  stability  margins  between  the  curves  a  and  /3  while  LTI  controllers 
are  useless,  it  may  be  the  case  that  suitable  nonlinear  control  provides  these  margins. 
Thus  unlike  the  case  of  completely  unstructured  modeling  uncertainty  (see  Khargonekar 
et  al.  [1987]  and  Pool  la  and  Ting  [1987])  the  possible  benefit  of  nonlinear  compensation 
cannot  as  yet  be  ruled  out  for  block-structured  problems.  This  interesting  open  issue  is 
currently  under  investigation. 

Kxamplc  f5.5).  Consider  again  an  additive  structured  family«of  plant  models  as  in  (2.3): 


V.v? 


0  A, 


A.  II  <51. 


It  follows  from  Theorem  (5.2)  that  the  maximal  stability  margin  R  for  the  family  above 
is  infinite.  Suppose  the  objective  is  to  construct  a  controller  K(s)  that  yields  a  stability 
margin  of  10  using  the  methods  of  Section  5.2.  Some  straightforward  calculations  yield 
the  following: 

Step  III:  y  =  10  .  R,  =  I2  and  R,  =  diag(4 , 0.25). 

Step  IV:  The  interpolation  problem  involved  is  to  find  an  outer  function  d(s)  such  that 
d(2)  =  4  and  d(3)  =  .25.  The  associated  Nevannlina-Pick  matrix  is 


N  = 

m 


4 m  H-  (.25) ' 


4  m  +(.25)m 


(.25)' 


and  it  can  be  verified  that  the  smallest  integer  m  for  which  Nm^0  is  m  =  7.  Thus,  one  can 
determine  a  14tb-order  outer  interpolant  via  (Krein  and  Nudelman  [1977])  to  be 


d(s)  = 


.6s“  +  322.8s  +  326.2  » 

130.2s2  +  134.8s  +  8.6 


Steps  V  and  VI:  The  desired  controller  can  now  be  obtained  as  in  Section  5.2  through 
tedious  buT  otherwise  straightforward  calculation. 
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CHAPTER  6 


SUBOPTIMAL  SOLUTIONS  FOR  MULTIDISK  PROBLEMS 


Recently,  the  single  objective  H^-optimization  problem  introduced  by  Zames  [1981] 
as  an  alternative  to  the  classical  Wiener-Hopf  approach  to  feedback  synthesis  has 
received  much  attention.  Several  interesting  and  important  control  problems,  however, 
involve  the  solution  of  a  multiobjective  H^-optimization  problem,  also  referred  to  as  a 


multidisk  problem  (see  Francis  and  Doyle  [1987]).  For  example, the  problems  of 
Robust  stabilization  with  optimal  nominal  disturbance  rejection 
Robust  simultaneous  stabilization 


(6.1) 

(6.2) 


Optimal  nominal  disturbance  rejection  with  robust  stability  around  a  failure  (6.3) 


operating  point 


all  reduce  to  (perhaps  nonlinear)  multidisk  problems. 

In  complete  generality  these  multidisk  problems  are  very  difficult  and  no  analytic 
solution  is  known.  This  chapter  shall  focus  instead  on  linear  two-disk  problems  (of  which 
(6.1)  is  an  example).  An  algorithm  will  be  presented  which  allows  computation  of  an 
upper  bound  for  linear  two-disk  problems  and  also  a  (suboptimal)  controller  which 
achieves  this  bound.  This  algorithm  readily  generalizes  to  incorporate  multidisk 
problems. 

This  chapter  is  organized  as  follows.  Section  6.1  contains  some  motivating  examples 
and  provides  some  background  information  about  two-disk  problems.  Section  6.2 
contains  ’he  main  results  of  this  chapter:  a  systematic  algorithm  for  computing  an 
aosoiute  upper  bound  for  a  linear  two-disk  problem  and  a  (subopumal)  controller  which 


s 


achieves  this  bound.  These  results  will  be  demonstrated  in  Section  6.3  through  a  few 
numerical  examples. 

6.1.  Background  Information  on  Multi-Disk  Problems 

In  Chapter  2  two  fundamental  problems  in  frequency  domain  control  were 
introduced:  the  robust  stabilization  problem  and  the  uniformly  or  H^-optimal  control 
problem.  Both  of  these  problems  involve  a  single  performance  criterion  and,  in  the 
context  of  LTI  controllers,  both  problems  reduce  to  the  same  mathematical  problem,  a 
one-disk  problem.  In  practice,  however,  a  control  system  designer  must  consider 
additional  performance  criteria,  i.e..  one  may  wish  to  simultaneously  consider  both  the 
robust  stability  properties  and  the  disturbance  rejection  capabilities  of  a  particular  control 
design.  For  example,  consider  the  system  shown  in  Figure  6.4  where  P  belongs  to  a  family 
of  plants  Fg  as  described  by  (2.2). 


FP-9290 


Figure  6.1.  Feedback  system  for  the  robust  stabilization  and 
uniformly  optimal  control  problem. 
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The  objective  is  to  design  a  controller  K  which  optimally  rejects  disturbances  d  for 
the  nominal  plant  P0  subject  to  the  constraint  that  K  robustly  stabilize  the  family  of 
plants  F5.  This  problem  is  equivalent  to  solving  the  constrained  optimization  problem 

min  II  XDW,  +  DQDW,  II  ^  subject  to  II  W,YD  -  W,NQD  II  ^  ^  1  , 

Q€RH„, 

which  represents  an  instance  of  the  two-disk  problem. 

In  complete  generality  multidisk  problems  are  very  difficult  and  no  analytic  solution 
is  known.  For  certain  special  classes  of  multidisk  problems,  a  partial  solution  has  been 
obtained  by  Ball  and  Helton  [1979].  Also,  Kwakernaak  [1982]  proposed  an  approximate 
solution  of  (6.1)  by  solving  instead  a  closely  related  single  objective  H^-optimization 
problem  ("mixed  sensitivity").  Finally,  O’Young  and  Francis  [1985]  explicitly  solve  a 
special  case  of  (6.1 )  under  some  reasonable  hypothesis. 

This  chapter  focuses  on  linear  multidisk  problems  (of  which  (6.1)  is  an  example). 
For  purposes  of  clarity  the  following  exposition  deals  only  with  two-disk  problems;  all 
results  easily  generalize  to  incorporate  more  than  two  competing  optimization  objectives. 
In  general  terms,  a  (linear)  two-disk  optimization  problem  may  be  formulated  as 
follows:  Given  matrix  Hw-  functions  A  ,  B  ,  C  ,  A  ,  B  ,  and  C  of  compatible  dimensions, 
determine  an  HM  function  Q  that  "jointly"  minimizes  II  A  +  BQC  B  and  H  A  +  BQC  I  M. 
Examples  of  this  joint  minimization  include 

min  II  A  4-  BQC  II  „  subject  to  II  A  +  BQC  B  „  <  y 
0 

where  y  is  a  prescribed  number,  or 

min  (  II  A  +  BQC  «  ^  +  9A  +  BQC  I  ^) 

Q 


v- 


(6.5) 


or,  more  generally, 


min  </>(  II  A  +  BQC  II  „  ,  IIA  +  BQCII^)  (66) 

0 

subject  to  \lK  II  A  +  BQC  II  „  ,  II  A  +  BQC  II  M)  <  1 
where  <b  and  \b  are  some  prescribed  functions. 

6.2.  The  Algorithm 

In  this  section  an  algorithm  is  presented  which  allows  computation  of  an  upper 
bound  for  linear  two-disk  problems  and  also  a  (suboptimal)  controller  that  achieves  this 
bound  (see  also  Ting  and  Poolla  [1987]).  This  algorithm  readily  generalizes  to  incorporate 
multidisk  problems.  Hara  and  Katori  [1986]  have  previously  investigated  the 
constrained  H^-optimization  problem  (which  is  an  instance  of  the  linear  two-disk 
problem).  For  single-input  single-output  systems  they  derive  upper  and  lower  bounds  on 
achievable  performance  bounds.  However,  their  lower  bounds  are  too  conservative 
because  they  do  not  incorporate  analyticity  requirements,  and  their  upper  bounds  are 
incomputable  because  they  involve  a  graph  metric  calculation  (see  Vidyasagar  [1984]). 
This  section  describes  a  computable  upper  bound  for  multiinput  multioutput  two-disk 
problems  which  requires  only  the  solution  of  several  one-disk  problems. 

To  begin,  consider  the  following 

Definition  (6.1):  Let  X  be  a  vector  space  and  let  f:X  -*  R+  and  G:X  -*  R+  be  nonnegative 
functionals  on  X.  Consider  the  problem  of  jointly  minimizing  (in  some  sense)  f  and  g 
over  X.  A  point  x  in  X  is  called  Pareto-optimal  if  there  does  not  exist  any  x  in  X  such 
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f(x)  ^  f(x*)  and  g(x)  <  g(x*) 

with  at  least  one  of  the  two  inequalities  being  si  net 

The  objective  of  this  definition  is  to  eliminate  clearly  suboptiinal  solutions  to  the 
problem  of  jointly  minimizing  f  and  g.  The  following  result  (see  DaCunha  and  Polak 
[l%7]  or  Vincent  and  Grantham  [1981])  characterizes  the  set  of  Pareto-optimal  solutions. 
Lemma  (6.1):  Let  X  be  a  convex  set.  and  let  f:X— »R+,  g:X  — >  R  +  be  convex 

s 

nonnegative  functionals.  Then,  x  in  X  is  Pareto-optimal  if  and  only  if  there  exists 
a  €  [0  ,  l]  such  that  x  minimizes  af(x)  +  (  1— u)  g(x)  over  X. 

□ 

Consider  now  the  linear  two-disk  problem  (6.4).  Observing  that  is  a  convex  set 
and  that  f(Q)=  II  A+BQC  II  and  g(Q)=  II  A  +  BQC  II  ^  are  convex  nonnegative 

functionals  on  HM  (via  the  triangle  inequality),  il  follows  that  Lemma  (6.1)  applies. 

Therefore,  the  problem  of  jointly  minimizing  f  and  g  can  be  solved  exactly  if  the  set  of 
all  Pareto-optimal  Q  can  be  computed.  For  instance,  consider  problem  (6.4).  For  each  a 
in  [0, 1  ]  suppose  it  is  possible  to  compute 

A  . 

arg  inf  [a  f(Q)  +  (  1  -  a)  g(Q)]  =Q  (a)  .  7) 

Q 

a  a 

Suppose  f(Q  (or))  and  g(Q  (a))  are  plotted  as  functions  of  a  as  in  Figure  6.2.  It  is  clear 
that  g  is  a  monotone  increasing  function  of  a  and  f  is  monotone  decreasing.  The 
constraint  g(Q)  ^  y  forces  one  to  consider  only  Pareto-optimal  Q  corresponding  to 

(»  -  [().<*„]  and  it  is  now  clear  that  an  optimal  solution  to  this  problem  is  Qopl  =  Q  (o-0). 


norm 


a  =0  alpha  a=l 

FP-9250 

Figure  6.2.  Plots  of  f(Q*(or))  and  g(Q*(or)). 

The  essential  difficulty  in  using  this  technique  to  solve  linear  two-disk  problems  is 
computing  Q  (a)  as  in  (6.7),  in  order  to  obtain  the  plots  in  Figure  6.2.  The  following  is  a 
step-by-step  description  of  an  algorithm  which  provides  a  (suboptimal)  approximation  of 
the  plots  of  Figure  6.2.  In  the  following  section  an  example  will  be  presented  which 
illustrates  this  algorithm. 

Consider  the  general  linear  two-disk  problem  (6.6). 

STEP  I:  Compute 

.  A  -  -  - 

Q  (0)  =arg  inf  II  A  +  BQC  II  „ 

QCH 

.  A  '  (6.8) 

Q  (  1 )  =arg  inf  II  A  +  BQC  II  M  . 

Q€H„, 

Remark  (6.1):  This  step  involves  only  the  solution  of  two  one-disk  problems,  which 
(Tom  considering  each  of  the  two  competing  objectives  individually.  Efficient 


computational  methods  exist  for  finding  solutions  to  these  problems  based  on  the  theory 
of  Hankel  operators  (see  Chu  and  Doyle  [1085],  Verma  [1985],  and  Glover  [1984]).  While 
the  solution  is  not  unique,  any  infemizing  choice  of  Q*  will  suffice. 


□ 


STEP  II:  For  a  in  [0.  l],  define 

A 

Q(a)=aQ(l)  +  (I  -a)Q(O)  (6.9) 

A  ..  A  - 

and  plot  f(or)  =  II  A  +  BQ(a)C  II  ^  and  g(a)  =  II  A  4-  BQ(a)C  II  ^  versus  a. 

Remark  (6.2):  f  and  g  serve  as  suboptimal  approximations  for  f  and  g.  it  readily  follows 
from  the  triangle  inequality  that  f  and  g  are  both  convex  functions  of  a. 


STEP  III:  Determine  the  set 


□ 


Q  =  [a  :  tiKf(a)  ,  g(a))  ^  I] 


Remark  (6.3):  In  this  step  one  may  use  these  approximations  to  determine  a  "feasible  set" 
U  of  a's  meeting  the  constraints. 


□ 

STEP  IV:  Compute  (graphically ) 

A  .. 

a0  =arg  inf  </>(f(a)  ,  g(a )) . 

at  Q 


Then,  a  (suboptimal)  upper  bound  solution  to  the  linear  2-disk  problem  (6.6)  is  given  by 

Qo  =  QK)- 

Remark  (6.4).  This  algorithm  yields  an  upper  bound  a0  because  the  choice  Q0  achieves 
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6.3  A  Numerical  Example 

In  this  section  the  algorithm  described  in  the  previous  section  will  be  demonstrated 
via  an  example.  Consider  the  family  of  plants  (as  in  Doyle  and  Stein  [1981]) 


I*g  =  {P0(I  +  AW,) :  II  A  II  <  1} 

and  all  plants  in  F5  are  constrained 


s+ 1 

(s+6)2 

s+2 

s2+2s+37 

s— 2  1 

where  P0(s)  = - ,  W,(s)  =  — 

s— 12  ‘  3 


to  have  exactly  one  unstable  pole.  The  design  objective  is  to  uniformly  optimally  control 
the  nominal  plant  P0  whilst  requiring  robust  stability.  More  specifically,  the  objective  is 
to  design  a  controller  K(s)  that 


i.  internally  stabilises  every  plant  in  Ffi,  and 

it.  infemizes  IIWj(I+P0K)  1  II  where  Wj(s)  =  — (-^—^-). 

30  s+1 

First  introduce  (as  in  Desoer  et  al.  [1980])  stably  coprime  factorizations  of  P0  as 
P0=N'D  1  with  XN+YD  =  1.  Using  standard  methods  one  obtains 


N  = 


s— 2 
s+6 


D  = 


s— 12 
s+6 


X  =  1.8  , 


Y  =  -0.8  . 


(6.10) 


Following  Theorem  (3.1 )  the  set  of  all  controllers  that  stabilize  P0  can  be  parameterized  as 

K(s)  =  (X  +  DQXY-NQ)-1  ,  Q  €  .  (6.11) 

From  the  well-known  results  of  Zames  and  Francis  [1983]  and  Doyle  and  Stein  [1981], 
the  design  problem  reduces  to 

inf  1  W,(I  +  P()Kf‘  I.  subject  to  1  W,P0K(I  +  P0K)~'  1.  <  1  .  ,,  p) 

K  =  M»bil.z;nj  P,  ' 

Substituting  ( 6. 1 1 )  above  yields 


ini'  II  W,  YD  -  W,NQD  II  „  subject  to  II  W,NX  +  W,NQD  0  „  <  1  . 
Qetd 


(6.13) 


which  is  an  instance  of  the  linear  two  disk  problem  (6.6)  with 


A  =  WjYD  ,  B  =  — W,N  .  C  =  D 


A  =  W,NX  ,  B  =  W,N.  C  =  D 


<t>=  II  A  +  BQC  II  „  ,  xls  =  II  A  +  BQC  II  M  . 

The  algorithm  described  in  Section  6.2  may  now  be  employed  to  obtain  a  suboptimal 
controller  for  this  feedback  synthesis  problem. 

STEP  I:  Using  standard  methods  the  following  solutions  to  the  one-disk  problems 
described  in  (6.8)  may  be  obtained 


LJ 


«  s— 21.3  »,  .  s+6 

Q(0)  =  0.2547  -  .  0  C 1 )  =  2.933  - 

s+1  s+12 


STEP  II:  One  may  now  plot  f(a)  and  g(o0  versus  a.  Numerically  this  results  in  the 
following  plot  of  Figure  6.3. 


m 
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STEP  III:  It  is  now  clear  that  the  set  Q  is 


Q  =  [  0 . 0.265  ] . 

STEP  IV:  a0  =  arg  inf  II  A  +  BQC  II  w  =  0.265  and 
«€  0 


Q0  =  .1872 


s— 21.3 
s+1 


+  0.7772 


s+6 
s+1 2 


Finally,  one  obtains  II  A  +  BQ0C  II  „  =  0.7998.  The  corresponding  controller  can  be 
determined  via  (6. 1 1 )  to  be 


1.8  + 

s— 12 

Qo 

-.8- 

s— 2 

Qo 

s+6 

s+6  , 

Remark  (6.5):  It  is  interesting  to  compare  the  upper  bound  for  the  linear  two-disk 
problem  with  the  results  achievable  when  considering  each  of  the  two  objectives 
individually.  For  example,  if  one  were  concerned  only  with  the  problem  of  uniformly 
optimal  control  of  P0  without  any  regard  for  the  robust  stability  of  IF5 ,  then  one  could 
reduce  the  H^-norm  of  the  weighted  sensitivity  function  to  0.1244.  Conversely,  if  one 

oo 

focuses  solely  upon  finding  the  optimal  robustly  stabilizing  controller,  then  the  H  -norm 
of  the  resulting  weighted  sensitivity  function  rises  to  1.043.  One  of  the  main  features  of 
the  curve  shown  in  Figure  6.3,  aside  from  providing  an  upper  bound  for  linear  two-disk 
problems,  is  to  graphically  display  the  design  tradeoffs  inherent  in  multi-objective 
problems. 


□ 


CHAPTER  7 


CONCLUDING  REMARKS 

This  dissertation  has  detailed  the  results  of  an  investigation  into  some  fundamental 
issues  in  control  system  design.  The  results  obtained  here  are  vital  in  answering  some 
basic  questions  regarding  what  can  and  cannot  be  accomplished  using  feedback  control. 
This  chapter  provides  a  brief  summary  of  the  main  results  and  contains  some  directions 
for  possible  future  research. 

In  Chapter  4  it  was  demonstrated  that  as  far  as  the  problem  of  robustly  stabilizing  a 
dynamic  family  of  plants  is  concerned,  the  best  controllers  are  linear  and  time-invariant. 
In  particular,  this  implies  that  adaptive  control  laws  (note:  these  generally  satisfy  the 
technical  hypothesis  of  Theorem  (4.5))  offer  no  advantage  over  LTI  controllers  in  the 
context  of  robustly  stabilizing  a  purely  dynamic  family  of  plants.  Due  to  the 
proliferation  of  available  controller  synthesis  methods,  results  of  this  type  are  invaluable 
in  ruling  out  various  synthesis  methods  for  certain  classes  of  control  problems. 

Adaptive  control,  however,  intends  primarily  to  deal  with  parametric  modeling 
uncertainty.  In  this  context,  an  interesting  result  due  to  Martensson  [1985]  demonstrates 
that  any  parametric  family  of  plant  models  with  compact  parameter  variation  can  be 
robustly  stabilized  by  using  some  nonlinear  time-varying  controller.  Martensson's 
treatment,  however,  deals  with  internal  (state)  stability  and  cannot  readily  be 
incorporated  into  the  input-output  framework  described  in  this  dissertation.  It  is 
nevertheless  apparent  that  for  robust  stabilization  problems  involving  purely  parametric 
uncertainty,  nonlinear  feedback  is  far  superior  to  linear  time-invariant  feedback. 
Blending  this  result  with  those  of  Chapier  4  one  is  led  to  conjecture  the  following  design 
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principle: 

The  precision  to  which  a  planl  model  is  known  determines  what  type  of  controller  is 

best  suited  for  controlling  the  plant. 

More  precisely,  it  appears  that  linear  time-invariant  controllers  are  optimal  for 
controlling  families  of  plants  featuring  primarily  dynamic  uncertainty.  Conversely,  it 
appears  that  families  of  plants  featuring  primarily  parametric  uncertainty  are  best 
controlled  by  nonlinear  time-varying  (or  adaptive)  controllers.  Intuitively,  this  follows 
from  the  belief  that  families  of  plants  featuring  primarily  parametric  uncertainty  could 
possibly  be  more  accurately  identified  through  further  input-output  experiments. 
However,  for  families  of  plants  featuring  primarily  dynamic  uncertainty,  further 
information  about  the  "true"  plant  model  cannot  be  obtained. 

A  representation  of  modeling  uncertainty  that  combines  both  dynamic  and 
parametric  uncertainty  is  the  motion  of  mixed  uncertainty  as  described  in  Chapter  2. 
This  representation  provides  a  setting  to  attempt  to  quantify  the  ideas  alluded  to  in  this 
design  principle.  In  other  words,  one  would  like  to  determine  how  much  dynamic 
uncertainty  can  be  exhibited  by  a  family  of  plant  models  before  nonlinear  time- varying 
controllers  no  longer  provide  significantly  better  performance  than  linear  time-invariant 
controllers.  Similarly,  one  would  also  like  to  determine  how  much  parametric 
uncertainty  can  be  contained  in  a  family  of  plant  models  before  nonlinear  time-varying 
controllers  become  far  superior  to  linear  time-invariant  controllers.  Recently,  (see,  for 
example.  Rohrs  et  al.  (l°82l)  there  has  been  interest  in  examining  how  adaptive  control 
laws  perform  with  unmodeled  dynamics.  One  possible  representation  of  this  modeling 
er-or  is  described  by  a  family  of  plants  featuring  mixed  uncertainty  as  given  in  (2.4),  i.e., 
1;5  represents  a  combination  of  parametric  and  dynamic  uncertainty.  An  interesting  open 
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question  then,  is  determining  an  exact  numerical  characterization  of  the  largest  radius 
8max  of  unmodeled  dynamics  that  can  possibly  be  stabilized  using  any  NLTV  controller. 

A  closely  related  problem  involves  the  construction  of  an  explicit  adaptive  or 
switching  control  scheme  for  a  family  of  plants  featuring  mixed  uncertainty.  More 
precisely,  consider  a  single  plant  P0  contained  in  a  family  of  plants  Fs  as  described  in 
(2.4).  Suppose  that  for  each  fixed  c*€  Q  there  exists  a  linear  time-invariant  controller 
that  robustly  stabilizes  the  resulting  family  of  plants.  The  key  question  then  is:  Does 
there  exist  some  method  of  gain-scheduling  or  switching  between  the  individual 
controllers  Kq  which  ensures  that  the  resulting  closed-loop  system  is  stable?  Along  these 
lines,  Fu  and  Barmish  [1986]  have  recently  applied  the  idea  of  switching  control  to  the 
problem  of  robustly  stabilizing  a  family  of  plants  with  compact  parameter  variations. 
They  exploit  the  fact  that  any  stabilizing  controller  for  a  given  nominal  plant  model 
contains  some  amount  of  robustness  to  develop  a  switching  law  between  a  finite  number 
of  controllers  which  guarantees  stability  for  the  true  plant  model.  Their  results  are  not 
directly  applicable  to  the  frequency  domain  robust  stabilization  problem  described  above 
because  they  deal  with  exponential  state  stability  as  opposed  to  input-output  stability. 
Such  problems  are  of  interest  to  the  adaptive  control  community  because  they  involve  a 
method  for  adaptive  stabilization  without  restrictive  assumptions  such  as  minimum-phase 
properties. 

In  Chapter  5  a  new  iterative  procedure  was  introduced  for  synthesizing  robustly 
stabilizing  controllers  for  families  of  plants  featuring  block-structured  uncertainties. 
This  method  is  referred  to  as  residue  iteration  and  provides  a  computationally  attractive 
alternative  to  the  established  A-Q  iteration  methods.  In  addition,  this  approach  was 
helpful  in  providing  valuable  insight  into  the  problem  and  triggered  the  development  of 


some  sufficient  conditions  on  the  nominal  plant  to  ensure  infinite  stability  margins.  To 
facilitate  the  computations  involved  in  performing  residue  iteration, a  new  method  was 
introduced  for  finding  low-order  outer  interpolating  functions.  This  new  procedure  is  also 
less  computationally  intensive  and  provides  lower-order  outer  functions  than  previously 
established  methods. 

Certain  numerical  aspects  of  the  residue  iteration  procedure  deserve  further 
attention.  Although  residue  iteration  is  numerically  attractive  as  compared  to  A-Q 
iteration,  both  procedures  share  a  common  difficulty  of  requiring  the  solution  of  a  non- 
convex  optimization  problem.  Recall  that  (5.1 )  is  not  jointly  convex  in  A  and  Q  and  (5.7) 
is  not  concave  in  the  residues  R_.  Thus,  neither  iterative  procedure  can  guarantee 
convergence  to  the  optimal  robustly  stabilizing  controller.  It  may  be  possible  to  show, 
however,  that  the  functional  a/j(M(R)W(R))  js  quasi-concave  in  the  residues  Rr  This 
would  be  desirable  because  quasi-concave  functions  have  at  most  one  local  maximum. 
Thus,  upon  arrival  at  a  local  maximum  in  (5.7)  it  is  only  necessary  to  check  the  boundary 
behavior  of  A.  (M(R)W(R))  as  R  approaches  oo  to  determine  the  global  maximum. 

The  development  of  residue  iteration  raises  several  related  issues  regarding  this 
robust  stabilization  problem.  For  example,  one  would  like  to  investigate  whether  it  is 
possible  to  extend  this  controller  synthesis  procedure  to  arbitrary  (non-square)  families 
of  plants  featuring  block-structured  uncertainty.  In  addition,  the  utility  of  nonlinear 
time-varying  controllers  for  robustly  stabilizing  families  of  plants  featuring  block- 
structured  uncertainty  is  a  completely  open  question. 

Chapter  6  detailed  a  step-by-step  algorithm  which  allows  for  computation  of  an 
upper  bourd  for  multiobjective  H^-optimization  problems.  These  results  are  obtained 
graphically  and  require  only  the  solution  of  several  single-objective  H^-optimization 


problems.  The  graphs  obtained  are  useful  in  demonstrating  the  design  tradeoffs  inherent 
in  multiobjective  problems. 


APPENDIX 


PROOF  OF  THEOREM  (4.3) 

Throughout  this  appendix,  M  is  a  (fixed),  m-input,  p-output  NLTV  operator  of  the 
form  M  =  A  +  BQC,  where  A,  B,  and  C  are  causal,  stable,  linear  time-invariant 
operators  with  C  being  rxr,  full  rank  and  Q  is  some  (fixed)  causal  incrementally  stable 
NLTV  operator.  Introduce  inner-outer  and  outer-inner  factorizations  of  B  and  C, 
respectively  (see  Averson  [1975])  as 

B  =  B;B0  .  C  =  C0C, 

Further,  assume  that 

BB*AC'lH  >  1  +  6,  ( A.  1 ) 

where  the  subscript  H  denotes  the  Hankel  norm  (see  Chapter  3). 

In  order  to  establish  the  proof  of  Theorem  (4.3)  the  following  intermediate  results 
shall  be  required. 

Lemma  (A.l).  Let  T  be  any  fixed  integer.  Then,  there  exists  an  input  sequence  w  in  h,m 
(dependent  on  T)  of  duration  K-T  and  such  that 

1  li  PKMw  I 

—  <  Iwl  <1,  P(w)  =  0,  -  >l  +  €.  (A. 2) 

2  lPKwl 

Proof.  First  note  that  from  Theorem  (4.2), 

inf  II  A  +  BQ,C  II  bh  =  II  B.'aC,’  II  h  >  1  +  c  .  (  x  3) 

:  causal, stableNLTV  3 

To  verify  (A. 2 ),  define  a  causal,  stable  NLTV  ooerator  Q,  :=  zTQz  t  where  z  1  is  the  unit 
delay  iperator.  It  now  follows  from  (A. 3)  that  there  exist  some  A  in  bh,  .uch  that 


86 


II  (A  +  BQ,C)\  II  /  II  K  II  >  1  +  e.  Define  v  =  z  K.  Then. 


II  Mv  II  II  (A  +  BQOz  T\  I!  II  zT(  A  +  BQC)z  T\  II 


II  (A  +  BQ,C)\  II 


>  I  +  €. 


Since  v  is  in  bh,m.  it  follows  that  there  exist  integers  T,  and  T,  such  that 


-  <  IIPNv  II  <  !  .  ail  N  ^  T, 

2 


II  Pn(A  +  BQC)v  II 


IIPNV 


>  1  +  e  .  all  N  >  T, 


Let  K  =  maxlTj.Til  and  define  w  =  PK(v).  It  is  now  easy  to  see  that  w  satisfies  (A.2), 
proving  the  claim. 


Proposition  (A.l).  Suppose  there  exists  some  input  sequence  w1  of  duration  Kj  such  that 


IPKM(w\v....)l 

o  Pk-  (w1. II 


>  1  +  €. 


Then,  there  exists  an  extension  w‘  of  duration  K2  —  K,  such  that 


—  <  II  w'  II  O  . 


I!  P..  M(w,w~,*, *,...)  1 


II  p  (w  ,w‘, *,*....)  II 


>  1  +  e. 


(A. 4) 


Proof.  Given  any  sequence  v  in  h,™  and  integers  K,.  K,  with  K,  <  K.-,.  define  the  finite 


sequence 


Let  <5  >  0  be  any  (fixed)  real  number.  Since  w  is  in  h,  and  since  M  is  stable,  it  follows 


rfoTrtAfa*- "Ur* • T  •  M ^ •  -  ’  •  K '  ’  ’’’ 


that  there  exists  an  integer  N  (without  loss  of  generality  N>  K.)  such  that 


U  -  PN  )M(  w  .0.0.  ■■■}!!  <  S. 


Let  A  =  C(  w  ,0.0,  ■  ■  ■  )  and  define  a  nonlinear  operator  Q  by 


()( u )  =  Q(  u  +  K )  —  Q(  A. )  . 


(A.5) 


It  is  easy  to  veril  y  that  (A.5)  defines  a  causal  operator.  Further,  since  Q  is  incrementally 
stable,  it  follows  that 

I!  Q(  u  +  A )  —  Q(  A )  II  ^ 

II  Q  II  =  sup  -  ^  II  Q  II  inc  <  oo  . 

Hud 

Thus,  Q  stable.  Define  M  =  A  +  BQC.  With  T  =  N  in  Lemma  (A.l)  it  follows  that 
there  exists  an  input  sequence  v  of  duration  K,  such  that 


-  ^  II  v  II  <  1  .  II  PNv  II  =  0  . 


II  P..  Mv 


II  Pv  v 


>  1  +  e. 


.cl  u"  =  [vl,  .  ,  and  notice  that  —  -Ss  II  w“  II  ^  1.  Notice  now  that 
N1'N2'  -) 


PK  M(w‘,w“,*, *,...)  II  ‘  =  II  PNM(w‘, II  '  +  II  (1  —  Pn)M(w\w‘, 


>  ( 1  +  e)“  H  w1  II  +  II  (PK  —  PN)  M(w\w",*,*,... )  II 


h(pKj-  PnInKw'.w2, 


^  II  p..  M(w’,w:, -  PK  M(  w1 ,0,0....)  a 


-  II  (IV  —  PvJM(  w  ,0,0,...)  I 
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>  II  Pv  Mv  H  -  8  ^  ( 1  +  e)  II  v  II  -  S  =  ( 1  +  e)  II  wx  II  -  8  . 


Combining  (A.6)  and  (A. 7)  yields 


1  Pv  M(w',w2, *,*,...)  II  2  >  ( 1  +  «)2(  II  w1  II  2  +  II  w2  II  2)  -  2S(  I  +  e)  II  w2  II . 


Since  8  can  be  chosen  to  be  arbitrarily  small  and  since  II  w2  II  ^  1,  it  follows  that  there 
exist  some  w"  such  that 

IIPkM(w\w2, *,*,...)  II  >  (I  +  £)!Pk(w1,w2iV,...)Ii 
K2  K2 


proving  ( A.4). 


One  is  now  in  a  position  to  prove  the  following 
Theorem  (4.3).  There  exists  an  input  sequence  u  in  h2™  such  that 


lim  II  PNu  II  =  oo 


(4.8) 


II  PNMu  II 
lim  - 

N-oo  II  PNU  II 


>  1  +  e. 


(4.9) 


Proof.  With  T  =  0  in  Lemma  ( A.l )  it  follows  that  there  exists  some  input  sequence  w1  of 
duration  K,  such  that 


-  <  II  w1  II  <  I  , 
o 


II  PK,M(w1. *,*,...)  II 
II  P.  ( w‘ ,*,*,...)  II 


>  1  +  €. 


(A.8) 


Now,  from  Proposition  (A.l),  one  may  conclude  that  there  exists  an  extension  w  of 


duration  K,  —  K,  such  that 
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II  PK  M(w\w2 
i<  lw!l  <1,  - ’ 


Pv  (w'.w2, 


>1  +  6. 


(A.9) 


Continuing  in  this  fashion,  it  is  possible  to  determine  using  Proposition  (A.l)  that  there 
exists  a  sequence  u  =  (w\w2,w3,...)  where  w1  is  of  duration  K;  --  Kj_, ,  with  the  following 
properties: 


-  <  II  w!  II  <  1  , 
o 


I  P^Mull 


H  PK  u  II 


>1  +  6. 


for  all  i . 


(A.10) 


Now  let  N  be  any  integer,  say  K;  <  N  ^  Ki+I.  Then. 


PnMu  II  11  pk  Mu  11 
- ^  — : - >  (1  +  €) 


PNu  II 


IIP.,  u  II 

Nl»l 


II  PK  u  B 


IIP,,  ull 

ki+i 


^  ( 1  +  e) 


•A 


i+2 


Consequently,  lim 


II  PNMu  II 


N— *oo  II  P|gU  II 


>1  +  6,  proving  (4.10),  and  clearly  u  is  unbounded. 


proving  (4.9).  This  completes  the  proof. 
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